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Introduction 


In  many  solid  tumors,  including  prostate  cancer,  hypoxia-inducible  factors  (HIF)  are  up-regulated 
compared  to  their  normal  tissue  counterparts.  These  HIF  molecules  are  transcription  factors,  and  up-regulate 
metabolic  and  angiogenic  proteins  in  the  cancer.  Recently,  a  set  of  proteins  known  as  prolyl  and  asparaginyl 
hydroxylases  (PHD  and  AHD  respectively)  have  been  shown  to  be  essential  in  the  regulation  of  HIF,  and  in 
some  cancers  have  been  transcriptionally  and  translationally  silenced.  We  therefore  proposed  a  study  that 
focuses  on  the  epigenetic  control  of  these  crucial  enzymes.  In  this  report,  we  present  data  demonstrating  our 
first  findings  of  PHD  expression  in  prostate  cancer  cell  lines  as  well  as  expanding  our  studies  to  relevant 
human  samples.  Furthermore,  we  begin  to  identify  specific  epigenetic  mechanisms  that  may  play  a  major  role 
in  the  transcriptional  and  translational  control  of  these  enzymes.  Last,  we  will  explain  our  future  direction  of  the 
project  after  the  award  period  has  expired. 

Body 

In  our  previous  reports,  we  had  begun  the  initial  characterization  of  the  PHD’s  and  AHD  in  an  array  of 
prostate  and  breast  cancer  cell  lines.  In  contrast  with  our  initial  hypothesis,  the  expression  levels  of  PHD1 , 
PHD2,  and  AHD  did  not  significantly  differ  amongst  the  cancer  cell  lines  compared  to  normal  controls.  On  the 
contrary,  the  mRNA  and  protein  level  of  PHD3  was  quite  variable  amongst  the  cell  lines  assayed.  Due  to  this 
observation,  our  initial  statement  of  work  was  changed  to  focus  on  this  specific  isozyme.  We  attempted  to 
examine  the  role  of  epigenetic  silencing  of  PHD3  in  these  cell  lines,  and  found  that  DNA  methylation  appeared 
to  be  a  key  regulator  of  the  expression  level  of  this  protein.  Furthermore,  it  was  found  that  PHD3  could  be 
reactivated  if  given  a  potent  DNA  de-methylating  agent  (5-aza-dC)  demonstrating  a  causal  relationship  to  the 
methylation  of  the  PHD3  promoter  and  gene  silencing.  Since  the  time  of  our  last  report,  we  have  confirmed  and 
extended  these  results  in  numerous  cancer  cell  lines  of  prostate,  breast,  and  melanoma  origins  (Figure  1).  The 
epigenetic  silencing  of  PHD3  was  not  limited  to  DNA  methylation,  as  it  was  observed  that  global  chromatin 
structure  was  also  altered  around  the  PHD3  promoter  (data  shown  in  previous  report).  While  these  results  are 
novel  and  exciting,  using  only  carcinoma  cell  lines  leaves  questions  about  clinical  relevance  of  our  findings.  To 
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address  this  issue,  we  have  finalized  studies  upon  prostate  tumor  samples  extracted  from  human  patients.  It 
was  discovered  that  on  average  the  clinical  tumors  demonstrate  significantly  lower  levels  of  PHD3  mRNA  than 
normal  prostate  tissue  (Figure  2).  Unfortunately,  the  limited  quantity  of  tumor  per  sample  we  were  able  to 
acquire  prevented  our  ability  to  confirm  these  results  at  the  protein  level.  In  contrast,  DNA  methylation  analysis 
was  performed  on  these  tissues  and  demonstrated  little  to  no  methylation  at  the  PHD3  promoter  regardless  of 
expression  level  of  PHD3  (data  not  shown).  Two  possibilities  exist  to  explain  this  finding:  1 )  The  region  within 
the  PHD3  promoter  that  was  analyzed  is  methylated  in  cancer  cell  lines  but  not  clinical  tumor  samples,  or  2) 
other  epigenetic  processes  may  be  involved  in  the  silencing  of  PHD3  in  vivo.  Both  of  these  hypotheses  provide 
enticing  endeavors  for  the  future  of  this  project. 

With  the  understanding  that  DNA  specific  epigenetic  marks  may  not  explain  the  entire  story  of  control  of 
PHD3,  we  attempted  to  examine  other  avenues  of  gene  expression  regulation.  When  examining  the  specific 
prostate  cancer  cell  line  DU  145,  it  was  observed  that  these  cells  had  attenuated  PHD3  expression  that  could 
be  reversed  by  hypoxia  but  exhibited  no  methylation  at  the  PHD3  promoter  (see  attached  manuscript  for  full 
data).  These  data  suggest  another  level  of  mRNA  regulation  beyond  DNA  methylation.  Using  the  miRNA 
prediction  software  TargetScan,  we  were  able  to  discover  4-5  potential  miRNA  binding  sites  in  the  PHD3  3’ 
untranslated  region  (UTR)  (Figure  3).  These  miRNA  sites  may  play  a  major  role  in  the  regulation  of  PHD3 
mRNA  levels  post-transcriptionally,  and  could  be  a  novel  regulatory  element  not  yet  described  for  this  gene.  To 
examine  if  these  miRNA  binding  sites  do  play  a  role  in  the  control  of  PHD3,  we  cloned  the  3’  UTR  of  PHD3  into 
a  vector  downstream  of  a  luciferase  gene.  By  doing  this,  we  allowed  the  stability  of  luciferase  to  be  controlled 
by  the  same  regulator  elements  that  may  impact  the  stability  of  PHD3.  Interestingly,  when  this  vector  was 
transfected  into  DU  145  prostate  cell  lines  it  was  observed  that  the  stability  of  luciferase  was  decreased  almost 
50%,  and  hypoxia  demonstrated  no  significant  change  in  this  reduction  (Figure  4).  While  these  data  are 
preliminary,  they  are  highly  suggestive  of  expression  regulatory  regions  within  the  3’  UTR  of  PHD3.  Future 
work  on  this  project  will  attempt  to  elucidate  the  exact  mechanism  by  which  PHD3  is  controlled  post- 
transcriptionally;  a  new  aspect  of  this  project  not  fully  appreciated  at  the  start. 
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Finally,  a  very  recent  discovery  has  shed  new  light  on  the  functional  role  of  PHD3.  Luo  et  al.  (Cell.  201 1 
May  27;145(5):732-44.)  have  shown  that  PHD3  may  interact  with  and  disrupt  the  function  of  pyruvate  kinase  in 
a  hydroxylase  independent  manner.  This  finding  implies  that  PHD3  may  have  a  direct  impact  on  cellular 
metabolism,  as  opposed  to  an  indirect  effect  through  the  hypoxia  inducible  factors  (HIFs).  To  query  if  this 
indeed  was  true,  we  examined  if  the  presence  of  PHD3  could  change  the  amount  of  oxidative  stress  (an 
indirect  measure  of  metabolism)  within  cancer  cells.  Using  an  adenovirus  to  over-express  PHD3,  we  infected 
the  pancreatic  cell  line  Mia-PaCa2  (this  cell  line  was  used  due  to  its  high  infectivity  with  adenovirus,  as  well  as 
undetectable  expression  of  endogenous  PHD3).  It  was  observed  that  the  over-expression  of  PHD3  in  these 
cell  lines  decreased  dihydroethidium  oxidation  by  approximately  half,  suggesting  a  50%  reduction  in  cellular 
reactive  oxygen  species  (Figure  5).  These  data  support  a  direct  metabolic  alteration  by  the  presence  of  PHD3, 
and  warrants  further  investigation  into  the  exact  mechanism  of  how  this  may  affect  tumor  cell  growth  and 
proliferation. 

In  summary,  we  have  thoroughly  characterized  the  epigenetic  processes  of  DNA  methylation  and 
chromatin  structure  at  the  PHD3  promoter  as  initially  outlined  in  our  statement  of  work.  While  certain  avenues 
of  this  project  did  not  produce  the  anticipated  results,  the  project  produced  significant  findings  in  the  regulation 
of  PHD3.  It  is  our  goal  to  continue  elucidating  novel  mechanisms  of  control  of  this  enzyme  in  the  hopes  of 
potential  new  therapies  for  patients  with  prostate  cancer. 

Key  Research  Accomplishments 

■  Completion  of  PHD3  expression  analysis  on  clinical  tumor  samples 

■  Completion  of  PHD3  promoter  methylation  analysis  on  clinical  tumor  samples 

■  Identification  of  putative  microRNA  binding  sites,  and  preliminary  regulatory  studies  using  the  3’UTR  of 
PHD3  on  a  luciferase  reporter 

■  Oxidative  stress  analysis  on  prostate  carcinoma  cell  line  in  the  presence  and  absence  of  PHD3,  with 
the  potential  for  understanding  metabolic  changes  due  to  PHD3 
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Conclusion 


In  summary,  this  project  has  surpassed  its  initial  goals  as  outlined  prior  to  the  funding  period.  Initially, 
our  study  was  limited  to  basic  expression  and  DNA  methylation  analysis,  but  we  have  expanded  our  abilities  to 
include  genetic  manipulations,  in-depth  oxygen  tension  studies,  miRNA  studies,  as  well  as  metabolic 
examinations.  This  work  has  produced  a  peer-reviewed  journal  article  in  a  highly  regarded  journal,  and  others 
are  in  preparation.  We  hope  to  continue  our  findings  on  PHD3  in  the  future,  and  anticipate  new  students  to 
take  on  the  challenges  of  elucidating  even  more  novel  findings  with  this  enzyme. 

The  PI  of  this  work  has  achieved  numerous  other  manuscripts,  research  awards,  and  has  presented  his 
work  at  a  variety  of  national  conferences  around  the  United  States.  Furthermore,  as  of  May,  201 1  and  the  end 
of  the  award  period,  the  PI  successfully  completed  his  PhD  studies  and  is  pursuing  a  1 -year  transition  post¬ 
doctoral  fellowship  with  Dr.  Domann  prior  to  pursuing  his  career  in  medicine/research. 


Case,  AJ  and  Domann,  FE 


Page  9 


Appendices  -  Figure  1 


PC-3  MCF7  MB-435  DU  145 


Figure  1.  Cell  lines  representative  of  PHD3  expressers  and  PHD3  non-expressers  were  either  not 
treated  (C),  treated  with  hypoxia  (H),  treated  with  5-Aza-dC  (Aza),  or  treated  with  5-Aza-dC  and 
then  subjected  to  1%  oxygen  (Aza  +  H).  Real  time  quantitative  PCR  was  performed  to  determine 
the  mRNA  expression  of  PHD3.  The  data  were  normalized  to  GAPDH  and  expressed  relative  to 
PHD3  mRNA  in  the  corresponding  untreated  controls.  Error  bars  =  SEM.  n  =  3. 
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Appendices  -  Figure  2 


A 


B 


Tumor  sample 


Figure  2.  PHD3  mRNA  expression  is  down-regulated  in  multiple  primary  human  prostate  cancer 
specimens.  A)  Total  mRNA  was  isolated  from  frozen  sections  of  primary  human  prostate  cancer 
specimens  with  Gleason  scores  ranging  from  7-9.  Quantitative  real-time  PCR  was  performed 
using  PHD3  specific  Taq Man  primer-probe.  Relative  PHD3  mRNA  expression  for  each  tumor 
sample  is  represented  as  (average  dCt  of  n  =  3  benign  prostate  tissue  samples)  -  (dCt  tumor 
sample).  Samples  were  normalized  to  GAPDH.  Dotted  lines  represent  +12  1  SD  for  benign  tissue 
PHD3  mRNA  expression.  B)  Box  plot  depicting  PHD3  mRNA  expression  from  samples  shown  in 
A.  p  value  is  based  on  ANOVA  between  3  benign  samples  and  10  malignant  samples. 
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Appendices  -  Figure  3 
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Figure  3.  TargetScan  analysis  of  3’UTR  of  PHD3.  When  analyzing  the  untranslated  region  of 
PHD3  four  major  putative  microRNA  binding  sites  are  elucidated. 
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Appendices  -  Figure  4 


DU145-  PHD3  3'UTR  Luciferase  Reporter 


Figure  4.  The  3’  UTR  of  PHD3  destabilizes  the  reporter  gene  luciferase.  After  successfully  cloning 
the  3’  UTR  of  PHD3,  the  segment  of  DNA  was  placed  downstream  of  a  luciferase  construct.  When 
transfected  into  DU  145  cells,  it  was  discovered  that  the  placing  luciferase  under  the  regulatory 
control  of  the  3’  UTR  demonstrated  an  approximate  50%  decrease  in  luciferase  stability. 
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Appendices  -  Figure  5 


MiaPaCa  -  DHE  analysis 
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Figure  5.  The  over-expression  of  PHD3  reduces  cellular  reactive  oxygen  species.  Mia-PaCa2 
cells  were  transfected  with  a  mammalian  expression  vector  of  PHD3,  and  then  exposed  to 
dihydroethidium  as  a  direct  measure  of  cellular  reactive  oxygen  species.  The  presence  of  PHD3 
reduced  dihydroethidium  staining  by  approximately  50%. 
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Abstract 

Background: The  prolyl-hydroxylase  domain  family  of  enzymes  (PHD1-3)  plays  an  important  role  in  the  cellular  response  to 
hypoxia  by  negatively  regulating  HIF-a  proteins.  Disruption  of  this  process  can  lead  to  up-regulation  of  factors  that  promote 
tumorigenesis.  We  observed  decreased  basal  expression  of  PHD3  in  prostate  cancer  tissue  and  tumor  cell  lines  representing 
diverse  tissues  of  origin.  Furthermore,  some  cancer  lines  displayed  a  failure  of  PFID3  mRNA  induction  when  introduced  to  a 
hypoxic  environment.  This  study  explores  the  mechanism  by  which  malignancies  neither  basally  express  PHD3  nor  induce 
PFID3  under  hypoxic  conditions. 

Methodology/Principal  Findings:  Using  bisulfite  sequencing  and  methylated  DNA  enrichment  procedures,  we  identified 
human  PHD3  promoter  hypermethylation  in  prostate,  breast,  melanoma  and  renal  carcinoma  cell  lines.  In  contrast,  non- 
transformed  human  prostate  and  breast  epithelial  cell  lines  contained  PHD3  CpG  islands  that  were  unmethylated  and 
responded  normally  to  hypoxia  by  upregulating  PHD3  mRNA.  Only  treatment  of  cells  lines  containing  PHD3  promoter 
hypermethylation  with  the  demethylating  drug  5-aza-2'-deoxycytidine  significantly  increased  the  expression  of  PFID3. 

Conclusions/Significance:  We  conclude  that  expression  of  PHD3  is  silenced  by  aberrant  CpG  methylation  of  the  PHD3 
promoter  in  a  subset  of  human  carcinoma  cell  lines  of  diverse  origin  and  that  this  aberrant  cytosine  methylation  status  is  the 
mechanism  by  which  these  cancer  cell  lines  fail  to  upregulate  PFID3  mRNA.  We  further  show  that  a  loss  of  PFID3  expression 
does  not  correlate  with  an  increase  in  HIF-1  oc  protein  levels  or  an  increase  in  the  transcriptional  activity  of  FHIF,  suggesting 
that  loss  of  PFID3  may  convey  a  selective  advantage  in  some  cancers  by  affecting  pathway(s)  other  than  FHIF. 
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Introduction 

The  cellular  response  to  reduced  oxygen  availability  (hypoxia)  is 
controlled  by  a  class  of  proteins  called  hypoxia-inducible  factors 
(HIF-a).  There  are  3  known  isoforms  of  HIF-la:  HIF-la,  HIF-2a 
and  HIF-3a.  HIF-la  and  ITIF-2a  are  transcription  factors.  HIF- 
3a  appears  to  lack  transcriptional  activity  and  may  play  a  role  in 
negative  regulation  of  the  HIF  pathway  [1],  Thus,  from  here  on, 
when  referring  to  HIF-a,  we  are  referring  to  only  HIF1  and  HIF2. 
Transcriptionally  active  HIF1  and  2  are  heterodimers  composed 
of  the  HIF-a  subunit  and  aryl  hydrocarbon  nuclear  translocator 
receptor  (ARNT/HIF-fl)HIF-la  activates  the  transcription  of 
EPO,  VEGF,  heme  oxygenase- 1  and  several  other  critical 


intracellular  responses  to  hypoxia  including  enzymes  of  the 
glycolytic  pathway  [2,3].  While  less  is  known  about  HIF-2a 
transcriptional  targets,  HIF-2a  appears  to  play  a  lesser  role  in  the 
glycolytic  response  with  more  emphasis  on  EPO  and  VEGF 
transcription  [4], 

HIF-a  mRNA  levels  are  generally  stable  in  cells.  It  is  not  until 
after  translation  that  HIF-a  is  tightly  controlled.  During  periods  of 
normal  physiological  oxygen  concentration,  HIF-a  subunits  are 
kept  at  low  levels  by  constant  proteolytic  degradation.  First,  a 
hydroxylation  reaction  is  catalyzed  by  a  family  of  prolyl 
hydroxylase  domain-containing  proteins  (PHD/EGLN/HPH) 
which  utilize  iron,  oxygen  and  2-oxoglutarate  as  co-factors  to 
enzymatically  catalyze  hydroxylation  on  the  oxygen-dependent 
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degradation  domain  (ODD)  of  the  HIFa  -subunit  [5].  Hydroxyl- 
ated  proline  residues  on  HIF-a  are  recognized  by  Von  Hippel- 
Lindau  (VHL)  protein,  an  E3  ubiquitin  ligase  that  ubiquitinates 
the  HIF-a  subunit,  targeting  it  to  the  proteosome  [6].  Under 
hypoxic  conditions,  HIF  prolyl  hydroxylase  activity  is  decreased 
and  HIF-la  protein  accumulates.  HIF-a  subunits  translocate  to 
the  nucleus  and  dimerize  with  the  constitutively  expressed  ARNT 
subunit  [7,8].  This  heterodimer  acts  to  turn  on  transcription  of 
genes  involved  in  oxygen  homeostasis  and  glucose  metabolism  [2] . 

Three  main  isoforms  of  HIF  prolyl-hydroxylase  domain 
containing  proteins,  PHD  1-3,  have  been  identified  [9].  These 
isoforms  have  been  reported  to  have  different  specificities  for  HIF- 
la  and  HIF-2a  [10],  and  also  differ  in  their  subcellular 
localization.  It  has  been  shown  that  PHD1  is  exclusively  present 
in  cytoplasm,  PHD2  is  mainly  located  in  the  nucleus  and  PHD3  is 
evenly  distributed  in  both  cytoplasm  and  nucleus  [1 1].  PHD2  and 
PHD3,  however,  are  considered  to  be  the  major  isoforms  that 
contribute  to  HIF-1  and  -2a  degradation  in  cells  [12,13].  In 
normoxia,  PHD2  is  the  primary  enzyme  that  hydroxylates  HIF-la 
[14],  whereas  PHD  3  has  been  reported  to  play  an  important  role 
in  HIF-2a  hydroxylation  and  also  in  retaining  cellular  hydroxyl- 
ation  capacity  in  a  hypoxic  environment  [10,15]. 

In  normal  cells,  PHD3  mRNA  and  protein  are  expressed  at  low 
levels  during  normoxia,  but  are  significantly  induced  upon 
exposure  to  hypoxia.  In  contrast,  PHD  3  expression  in  a  significant 
number  of  cancer  cell  types  has  been  shown  to  be  low  or  absent 
not  only  during  normoxia,  but  also  under  hypoxic  conditions 
[10,16].  To  date,  no  mechanism  has  been  discovered  to  explain 
this  defect  in  hypoxic  inducibility.  Interestingly,  Hatzimichael  et 
al.  have  recently  demonstrated  that  the  promoter  of  PHD  3  is 
methylated  in  certain  primary  B-cell  dyscrasias  [17].  We  had 
observed  a  decrease  in  PHD3  mRNA  expression  in  human  breast 
and  prostate  carcinoma  cell  lines,  with  an  absence  of  PHD3 
upregulation  in  response  to  hypoxia.  Therefore,  we  were 
interested  to  determine  whether  PHD  3  promoter  methylation 
was  responsible  for  this  aberrant  expression  pattern.  In  this  study, 
we  show  that  the  promoter  region  of  PHD3  is  mediylated  in 
representative  human  prostate  carcinoma,  melanoma,  renal 
carcinoma  and  breast  cancer  cell  lines.  Furthermore,  we  show 
that  neither  HIF-la  protein  levels  nor  hypoxic  response  through 
an  HRE-luciferase  reporter  vector  are  compromised  in  PHD  .3 
methylated  compared  to  non-methylated  cell  lines.  These  results 
indicate  that  PHD  3  promoter  methylation  is  utilized  by  malig¬ 
nancies  derived  from  diverse  human  cell  types.  Furthermore,  these 
data  suggest  that  loss  of  PHD  3  expression  may  not  affect  the 
transcriptional  response  through  the  HIF  pathway,  leaving  open 
the  possibility  that  PHD3  silencing  in  tumors  is  selected  through 
the  loss  of  specific  interactions  with  other  cellular  pathways. 

Methods 

Cell  culture 

Normal  human  prostate  epithelial  cells  (NPrEC)  were  pur¬ 
chased  from  Clonetics,  Lonza  Inc.  (Walkersville,  MD)  and  were 
grown  on  the  recommended  PrEGM  media  supplied  by  Clonetics, 
Lonza  Inc.  The  hTERT-HMEl  cells  were  cultured  in  mammary 
epithelial  basal  medium  MEGM  (Lonza  Inc.)  at  37°C  and  5% 
CO2  according  to  the  manufacturer’s  instructions  (Lonza  Inc.). 
The  DU  145,  22RV.1,  PC-3,  MDA-MB-435  (MB-435),  and 
MCF7  cell  lines  were  obtained  from  ATCC  (Manassas,  VA). 
MCF7,  DU  145  and  MDA-MB-435  cells  were  cultured  in  Eagle’s 
Minimum  Essential  Medium  (MEM)  supplemented  with  10%  fetal 
bovine  serum  (FBS).  PC-3  cells  were  cultured  in  FI 2  medium 
supplemented  with  10%  FBS,  2  mM  L-glutamine,  1  inM  Na 


Pyruvate,  and  supplemented  with  100  U/ml  pen/strep.  All  cell 
lines  were  routinely  maintained  at  37°C  in  a  humidified 
atmosphere  with  5%  C02.  Fresh  media  was  replaced  every  three 
days  while  routine  subculture  was  performed  by  washing  with  IX 
PBS  and  detaching  cells  with  TrypLE  Express. 

Semi-quantitative  RT-PCR 

Total  RNA  was  extracted  from  individual  cell  lines  using  RNeasy 
Mini  Kit  (Qiagen,  Valencia,  CA)  and  quantified  using  a  NanoDrop 
1 000.  To  assess  PHD3  and  GAPDH  expression,  500  ng  of  total  RNA 
was  used  for  reverse  transcription  using  a  OneStep  RT-PCR  Kit 
(Qiagen).  The  PHD3  forward  primer  is  S'-GGGCAAATAC- 
TACGTCAAGGAG-3'  and  the  reverse  primer  is  5'-AGTCTT- 
CAGTGAGGGCAGATTC-3'.  GAPDH  expression  was  assessed 
using  GAPDH-specific  primers.  PGR  conditions  for  PPID3  and 
GAPDH  were  the  same  except  that  28  cycles  of  PCR  were  performed 
for  PHD3  analysis  and  23  cycles  were  performed  for  GAPDH.  The 
parameters  used  were:  95°C  for  5  minutes  followed  by  the  stated 
number  of  cycles  of  94°C  for  1  minute;  56°C  for  1  minute,  and  72°C 
for  1  minute,  ending  with  a  final  extension  at  72°C  for  7  minutes.  The 
amplified  products  were  electrophoresed  on  a  1%  agarose  gel  and 
stained  with  ethidium  bromide  to  visualize  the  bands. 

Quantitative  real  time  RT-PCR 

Total  RNA  was  isolated  from  cells  using  Trizol,  followed  by 
DNAse  treatment  and  NaOAc  precipitation.  The  reverse  transcrip¬ 
tion  reaction  was  carried  out  with  High-Capacity  cDNA  Archive 
Kit  (Applied  Biosystems,  Foster  City,  CA).  PHD3  TaqMan  primer- 
probe  was  utilized  from  Applied  Biosystems  (Hs00222966_ml). 
The  quantitative  real-time  PCR  was  set  up  as  follows:  10  ngofRNA 
was  used  as  template  for  each  real-time  PCR  reaction  (10  pg 
reaction  volume);  primer  pairs  at  0.3  pM  for  GAPDH  with  Syber 
Green  Master  Mix  (Applied  Biosystems).  For  PHD3,  TaqMan 
universal  master  mix  was  used.  The  DNA  polymerase  was  activated 
by  heat  at  95°C  for  10  min  followed  by  40  cycles,  denaturing  at 
95°C  for  15  s,  annealing  and  elongating  at  60°C  for  1  min.  Data 
were  collected  with  ABI  PRISM  7000  sequence  detection  system. 
Data  were  analyzed  using  the  AACt  method. 

Western  Blot  analysis 

Cells  were  immediately  washed  with  ice  cold  phosphate- 
buffered  saline  (pH  7.4).  Cells  were  lysed  on  the  plate  in  200  pi 
RIPA  cell-lysis  buffer  (50  mM  Tris  pH  8.0,  150  mM  NaCl,  0.1% 
SDS,  0.5%  Na  Deoxycholate,  1%  TX-100)  plus  1  mM  NaF, 
10  mM  NaVC>4,  10  mM  PMSF,  and  1/100  protease  inhibitor 
cocktail  (Sigma),  immediately  boiled  for  2  minutes  and  then 
sonicated.  SDS-polyacrylamide  gels  (7%,  PHD3;  15%  HIF1, 
HIF2)  were  used  for  protein  electrophoresis.  Proteins  were 
electrotransferred  onto  nitrocellulose  membranes  and  treated  with 
anti  HIF-la  (Abeam,  Cambridge,  MA)  1:500  overnight  at  4°C. 
Anti  PHD3,  NB 100- 139  and  anti  HIF-2a  antibodies  (Novus 
Biologicals,  Littleton,  CO)  were  used  at  1:500  and  1:200 
respectively  overnight  at  4°C.  Equal  protein  loading  was 
confirmed  on  all  immunoblots  using  human  actin  antibody 
(Sigma,  St.  Louis,  MO)  at  a  dilution  1:2000.  Goat  anti-rabbit 
IgG  (BD  Transduction  Laboratories,  San  Diego,  CA)  was  used  as 
a  secondary  antibody  against  all  primary  antibodies.  Bands  were 
visualized  by  chemiluminescence  with  ECL  plus  reagent  (Pierce, 
Rockford,  IL)  on  a  Typhoon  FLA  7000. 

Sodium  bisulfite  sequencing 

Genomic  DNA  was  extracted  with  the  use  of  the  DNeasy  Tissue 
Kit  (Qiagen,  Valencia,  CA),  and  sodium  bisulfite  conversion  was 
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performed  with  the  use  of  the  EZ  DNA  Methylation  Kit  (Zymo 
Research  Corporation,  Orange,  CA).  A  pair  of  primers  was 
designed  to  amplify  the  PHD3  promoter  of  both  bisulfite  modified 
methylated  and  unmethylated  DNA  but  not  unmodified  DNA. 
Nested  PCR  amplification  on  converted  DNA  used  the  following 
primers:  outside  foiward:  5'-GTGTGGGATTTAGGTTTT- 
TAAG-3'  (SB1);  outside  reverse:5'-CCAAATCCAACCTCA- 
TAATATATC-3'  (SB2);  and  nested  inner  primers  (SB3)  and 
(SB4)  whose  sequences  and  locations  are  described  in  detail  below. 

The  resulting  PCR  products  were  gel-extracted  with  the  use  of 
the  Qiagen  Gel  Extraction  Kit,  or  gel  digestion  with  (3-agarase 
followed  by  EtOH  precipitation,  and  cloned  with  the  TOPO  TA 
Cloning  Kit  (Invitrogen).  Plasmid  DNA  was  extracted  with  the  use 
of  the  QiaPrep  Spin  Plasmid  Miniprep  Kit  (Qiagen).  Sequencing 
was  performed  by  the  sequencing  core  facility  maintained  by  the 
University  of  Iowa  and  results  were  tabulated  for  methylation 
status  of  each  of  the  58  CpGs  contained  in  the  amplicons  from 
each  cell  line. 

5-Aza-dC  treatment 

Cells  were  counted,  and  seeded  (day  0)  at  approximately 
750,000  cells/ 100  mm  dish.  Fresh  5-Aza-dC  (5  mM)  was  added  to 
the  dish  on  days  1,  3,  and  5  while  a  control  flask  was  left  untreated. 
On  day  5,  5-Aza-dC  treated  cells  were  split  into  two  60  mm  dishes 
in  media  supplemented  with  5  mM  5-Aza-dC.  On  day  6,  one  of 
each  of  the  60  mm  dishes  was  placed  in  a  hypoxia  chamber  and 
placed  under  1%02,  94%N2  and  5%  C02  at  37°C.  On  day  7,  all 
cells  were  harvested  with  500  pi  Trizol  for  RNA  extraction. 

Chromatin  Accessibility 

Chromatin  accessibility  experiments  were  conducted  as  previ¬ 
ously  described  by  Rose  et  al.  [18].  The  primers  CA1  and  CA2 
were  located  within  the  region  of  the  PHD3  gene  queried  for  DNA 
methylation  (primer  sequences  and  locations  described  in  detail 
below).  After  nuclei  extraction  and  a  5  minute  DNase  I  digestion 
the  DNA  was  extracted  and  real-time  PCR  was  conducted  on  an 
ABI  7000  Sequence  Detection  System.  The  accessibility  index  for 
each  amplicon  was  then  determined  by  the  following  formula 
(accessibility  index  =  2((Q  Dnase  treated)  “  (Q  Uncu,))).  GAPDH 
chromatin  accessibility  was  also  determined  as  a  positive  control 
for  a  constitutively  expressed  gene  to  control  for  equivalent  DNase 
digestion  between  the  cell  lines  examined. 

HRE-Luciferase  assay 

Cell  lines  ~85%  confluent  in  60  mm  dishes  were  transfected 
with  an  HRE-luciferase  reporter  vector  [19]  (2.5  pg)  and  Renilla 
luciferase  (1.5  pg)  according  to  Lipofectamine  2000  transfection 
reagent  protocol.  Transfection  media  was  removed  after  6  hours 
and  replaced  with  fresh  medium.  Cells  we  then  placed  under  94% 
N2,  5%  C02,  1%  02  gas  mixture  in  a  hypoxia  chamber,  or 
normoxia  for  24  hours  and  then  lysed  according  to  the  Dual 
luciferase  reporter  assay  system  (Promega,  Madison,  WI)  protocol. 
Luminescence  was  measured  3  times  per  sample  using  a  Tecan 
SpectraFluor  Plus  luminometer. 

Adenoviral  Transduction 

The  adenoviral  PF1D3  construct  was  a  generous  gift  from  Dr. 
Robert  Freeman  from  Rochester  University  (unpublished). 
Briefly,  the  human  PHD3  coding  sequence  was  engineered 
into  the  pDC315  vector  and  contains  an  N-terminal  FLAG  tag. 
PC3  cells  were  grown  to  ~85%  confluency  and  then 
transduced  with  20,  40  or  60  MOI  of  Ad-PHD3.  Approxi¬ 
mately  36  hours  following  transduction,  cells  were  lysed  with 


RIPA  buffer  and  western  blotted  according  to  the  procedures 
outlined  above. 

Clinical  Samples 

Clinical  prostate  tumor  samples  were  received  as  frozen  blocks 
in  OCT.  Sections  were  cut  and  ground  with  mortal  and  pestle. 
DNA  and  RNA  were  extracted  in  with  Qiagen  DNeasy  Tissue  kit 
and  Trizol  respectively. 

Methylated  DNA  Enrichment 

Genomic  DNA  was  harvested  from  cells  and  tissue  using  a 
Qiagen  DNeasy  Tissue  kit.  2  pg  of  DNA  in  120  pi  of  10  mM  Tris 
pH  8.0  was  sonicated  into  fragments  of  approximately  150  bp 
using  a  Covaris  S2.  Fragmentation  was  done  according  to  the 
Covaris  protocol.  1  pg  (60  pi)  of  sheared  input  gDNA  was  used 
according  to  the  protocol  supplied  by  the  MethylMiner  kit 
(Invitrogen).  Binding  reactions  between  beads  containing  methyl- 
CpG  binding  domains  and  sheared  genomic  DNA  were 
performed  at  4°C  overnight.  Bound  DNA  was  eluted  using 
progressively  increasing  NaCl  concentrations.  Eluates  were 
precipitated  using  NaOAc  and  EtOH  precipitation  and  resus¬ 
pended  in  60  pi  10  mM  Tris  pH  8.0.  Real  time  PCR  was 
performed  using  1  pi  resuspended  DNA,  SYBR  Green  master  mix 
and  100  nM  Fwd  Primer:  5'-GAGCTCCACGACCCGTTTC-3' 
and  Rev  Primer:  -5'-GCAGTGGTGGCTTCCCAT-3,  in  a  10  pi 
reaction  volume.  The  kit  was  validated  using  samples  from  human 
tumor  cell  lines  with  know  PHD3  CpG  island  methylation  status  as 
determined  by  bisulfite  sequencing  (Figure  SI). 

Statistical  Analysis 

Significant  differences  between  groups  of  data  were  determined 
using  a  t-test  for  all  bar  graphs  or  ANOVA  for  box  plot  powered 
by  SigmaPlot  1 1 .0  software,  n  =  3  was  used  for  each  data  set  unless 
otherwise  noted. 

Results 

PHD3  mRNA  is  aberrantly  silenced  in  human  melanoma, 
prostate  and  breast  carcinoma  cell  lines 

A  panel  of  human  carcinoma  cell  lines  was  screened  for  PHD  3 
mRNA  expression.  This  panel  consisted  of  three  prostate  cancer 
cell  lines  (DU  145,  22RV.  1  and  PC-3),  two  breast  cancer  cell  lines 
(MCF7,  HS578T),  one  melanoma  (MDA-MB-435),  a  non- 
transformed  prostate  epithelial  cell  line  (NPrEC)  and  a  non- 
transformed,  immortalized  breast  cell  line  (hTERT-HMEl).  We 
found  that  PHD3  mRNA  was  expressed  at  different  levels  varying 
from  abundant  to  almost  undetectable  levels  as  determined  by 
conventional  RT-PCR  (Figure  la).  The  prostate  cancer  cell  lines 
showed  decreased  PHD3  mRNA  expression  compared  to  the 
normal  prostate  epithelial  cells.  A  comparison  of  three  prostate 
cancer  cell  lines  showed  that  PHD3  is  expressed  in  DU  145  and  22 
RV.  1 ,  whereas  in  PC-3,  PHD3  mRNA  is  nearly  undetectable.  The 
melanoma  cell  line,  MDA-MB-435  was  also  found  to  have  very 
low  PHD  3  mRNA  expression.  Among  the  mammary  cell  lines, 
HS578T  had  much  lower  PHD3  mRNA  than  the  HME1 
mammary  epithelial  cells  and  MCF7  cells  expressed  far  more 
PHD  3  mRNA  than  its  normal  HME1  mammary  epithelial  cell 
counterpart.  To  confirm  and  extend  the  results  shown  in 
Figure  la,  quantitative  real  time  RT-PCR  analysis  of  PHD3 
mRNA  expression  were  conducted  and  the  results  are  shown  in 
Figure  lb.  Similarities  in  expression  were  found  with  both  the 
methods  used,  and  PHD  3  mRNA  expression  was  nearly 
undetectable  in  PC-3,  MDA-MB-435  and  HS578T  cell  lines. 
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Figure  1.  Differential  expression  of  PHD3  mRNA  in  normal  and  tumor  cell  lines  of  prostate,  breast  and  melanoma  origin.  Total  RNA 
was  extracted  from  the  indicated  cell  lines.  A)  Expression  of  mRNA  was  determined  by  semiquantitative  RT-PCR  using  specific  primers  designed  to 
amplify  PHD3  and  GAPDH.  The  RT-PCR  products  were  electrophoresed  using  1%  agarose  gel,  with  ethidium  bromide  to  show  bands  of  expected 
sizes.  GAPDH  was  used  as  the  loading  control.  B)  Quantitative  real-time  reverse  transcription-PCR  analysis  of  PHD3  was  performed  with  normalization 
to  18S  rRNA  gene  expression  and  relative  quantities  were  determined  by  the  DDCt  method.  The  data  shown  are  the  PHD3  gene  expression  of 
individual  cancer  cell  lines  relative  to  those  of  their  corresponding  normal  counterparts  (PrEC  for  prostate  and  HMEC-hTERT  for  mammary). 
doi:10.1371/journal.pone.0014617.g001 


Re-expression  of  PHD3  with  DNA  methyltransferase 
inhibitor  5-Aza-dC 

The  near  absence  of  PHD3  mRNA  expression  in  a  subset  of  cell 
lines  suggested  an  epigenetic  mechanism  might  be  responsible  for 
their  silencing.  Unlike  genetic  mutations  that  accumulate  in 
cancer,  epigenetic  modifications  are  reversible  [20].  We  hypoth¬ 
esized  that  if  DNA  methylation  of  the  PHDS  gene  was  responsible 
for  its  reduced  expression,  a  DNA  methyltransferase  inhibitor  such 
as  5-Aza-2'-deoxycytidine  (5-Aza-dC)  should  induce  its  expression. 
This  is  in  line  with  5-Aza-dC ’s  purported  ability  to  re-activate 
genes  previously  silenced  by  DNA  methylation  in  cancer  cells  [21]. 
We  chose  2  cell  lines  displaying  the  most  marked  decrease  in 
PHD3  expression,  MB-435  and  PC-3,  and  2  cell  lines  displaying 
moderate  to  high  basal  PHD3  expression,  MCF7  and  DU  145  for 
treatment  with  5-Aza-dC.  When  MB-435  and  PC-3  cells  were 
treated  with  5-Aza-dC,  there  were  significant  increases  in  the 
PHD3  mRNA  expression  compared  to  their  respective  untreated 
controls  (Figure  2).  Furthermore,  MB-435  cells  became  respon¬ 
sive  to  PHD3  mRNA  upregulation  by  hypoxia  to  a  significant 


degree  following  hypoxic  exposure.  Cells  that  already  expressed 
PHD  3  at  moderate  levels  did  not  respond  to  5-Aza-dC  by 
significantly  upregulating  PHD3.  Although  this  finding  suggests 
that  CpG  methylation  is  involved  in  silencing,  further  direct 
queries  of  epigenetic  alterations  at  this  locus  were  necessary  to 
more  deeply  address  this  question. 

The  PHD3  CpG  island  is  aberrantly  CpG  methylated  and 
displays  decreased  chromatin  accessibility  in  human 
carcinoma  cells 

T o  determine  whether  PHD3  gene  methylation  is  present  at  the 
CpG  island  in  cell  lines  that  have  reduced  PHD  3  mRNA 
expression  and  respond  to  5-Aza-dC  by  upregulating  PHD  .3 
mRNA,  we  utilized  sodium  bisulfite  sequencing  to  identify 
methylated  CpG  sites.  Figure  3A,  illustrates  the  CpG  island  in 
the  5*  -end  of  the  PHD3  gene.  The  58  CpG  sites  in  the  region 
analyzed  are  represented  by  vertical  lines.  A  putative  hypoxia 
response  element  (HRE)  in  the  promoter  region  of  the  gene  is  also 
indicated.  We  found  that  these  CpGs  were  highly  methylated  in 
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Figure  2.  Treatment  with  5-Aza-dC  triggers  re-expression  of  PHD3.  Cell  lines  representative  of  PHD3  expressers  and  PHD3  non-expressers 
were  either  not  treated  (C),  treated  with  hypoxia  (H),  treated  with  5-Aza-dC  (Aza),  or  treated  with  5-Aza-dC  and  then  subjected  to  1%  oxygen  (Aza  + 
H)  according  to  the  protocol  described  in  Materials  and  Methods.  Real  time  quantitative  PCR  was  performed  to  determine  the  mRNA  expression  of 
PHD3.  The  data  were  normalized  to  GAPDH  and  expressed  relative  to  PHD3  mRNA  in  the  corresponding  untreated  controls.  Error  bars  =  SEM.  n  =  3. 
doi:1 0.1 371  /journal.pone.001 461 7.g002 


the  PHD3  negative  cell  lines  PC-3,  MB-435,  and  HS578T 
(Figure  3B).  These  represent  examples  of  human  melanoma, 
prostate,  and  mammary  carcinoma  cells  respectively.  We  also 
noted  that  many  cell  lines  appeared  to  be  heterogeneous  with 
respect  to  PHD3  promoter  methylation  status.  Within  certain  cell 
lines,  some  clones  display  high  levels  of  methylation,  whereas  other 
have  very  few  to  no  methylated  CpGs.  This  small  population  of 
unmethylated  or  hemimethylated  cells  within  a  cell  line  may 
explain  our  ability  to  detect  very  low  levels  of  PHD  3  mRNA  in  cell 
lines  displaying  largely  methylated  PHD3  CpG  islands.  Further¬ 
more,  areas  of  CpG  methylation  in  some  methylation  positive  cell 
lines  overlap  with  a  putative  HRE  in  the  PHD3  promoter  region, 
which  could  hinder  the  ability  of  PHD  3  to  be  induced  upon 
hypoxic  stimuli.  In  contrast  to  the  positive  methylation  status  in 
PHD3  negative  cells,  the  CpGs  in  the  PHD3  CpG  island  were 
largely  unmethylated  in  the  PHD3  positive  cells,  NPrEC,  DU  145, 
22RV.1,  HME 1 ,  and  MCF7. 

DNA  methylation  is  typically  associated  with  other  alterations  to 
chromatin  structure  that  participate  in  cell-type  specific  gene 
expression  patterns.  Aberrant  cytosine  methylation  in  the  5'- 
regulatory  regions  of  genes  is  typically  associated  with  deacetylated 
histones,  and  thus  a  state  of  DNA  that  is  generally  inaccessible  to 
transcription  factors  and  other  enzymes  that  act  on  DNA,  such  as 
polymerase  II.  This  is  a  mechanism  of  gene  silencing  often 
exploited  by  cancer  cells  [22].  In  a  chromatin  accessibility  assay, 
we  found  the  promoter  region  of  PHD3  in  PC-3  cells  was  resistant 
to  cutting  by  DNase  I  when  compared  to  the  MCF7  PHD3 
promoter,  whereas  there  was  little  change  in  GAPDH  promoter 
accessibility  between  the  two  cell  lines  (Figure  S2).  This  evidence 
further  supports  the  hypothesis  that  PHD3  promoter  mediylation 
and  heterochromatin  formation  are  part  of  the  mechanism  for 
reduced  expression  of  PHD3  in  these  human  breast  and  prostate 
cancer  cell  lines. 


PHD3  expression  is  not  induced  upon  exposure  to 
hypoxia  in  cell  lines  containing  PHD3  promoter 
methylation 

Unlike  PHD1,  both  PHD2  and  PHD3  genes  contain  hypoxia 
response  elements,  and  can  be  induced  by  hypoxia  by  the  HIF-1 
and  HIF-2  transcription  factor  complex.  In  the  case  of  PHD3, 
mRNA  and  protein  expression  can  be  relatively  low  during 
normoxic  conditions,  with  marked  increases  upon  hypoxic  insult 
[10].  Therefore,  we  tested  PC-3,  DU  145,  MB-435  and  MCF7  cell 
lines  for  their  ability  to  upregulate  PHD3  following  24  hours  of 
hypoxia  (1%  O2)  (Figure  4A).  We  found  cell  lines  that  contained 
PHD3  promoter  methylation  (PC-3,  MB-435)  failed  to  appreciably 
upregulate  PHD3  mRNA  under  these  conditions.  However,  we 
did  note  a  very  small  upregulation  of  PHD3  mRNA  in  PC-3  cells. 
This  can  likely  be  attributed  to  the  heterogeneity  of  PHD3 
promoter  methylation  between  specific  clones  in  this  cell  line  (see 
figure  3B).  In  contrast,  PHD  3  mRNA  was  much  more  prone  to 
upregulation  in  the  unmethylated  cell  line  MCF7.  Upregulation  of 
PHD3  in  DU  145  cells  varied  by  experiment,  and  averaged  as  a 
non-significant  trend  toward  hypoxic  upregulation. 

Our  discovery  of  PHD3  promoter  methylation  in  melanoma, 
breast  and  prostate  cancer  cell  lines  prompted  us  to  ask  whether 
cell  lines  from  other  malignant  tissues  contained  methylation  at  the 
PHD3  locus.  Therefore,  we  further  performed  real-time  PCR  and 
bisulfite  sequencing  on  a  panel  of  3  human  renal  cell  carcinoma 
cell  lines  (A-498,  ACHN  and  769-P),  and  on  cDNA  prepared  from 
normal  kidney  tissue  (Figure  4B).  We  found  that  ACHN  and 
769-P  cells  express  nearly  undetectable  levels  of  PHD  3  mRNA, 
whereas  A-498  expresses  levels  comparable  to  normal  tissue.  Of 
these  cell  lines,  the  PHDS  positive  A-498  displayed  an  unmethy¬ 
lated  PHD  3  promoter  whereas  769-P  cells  displayed  an  aberrantly 
methylated  PHD  3  promoter.  Interestingly  though,  we  did  not 
detect  any  CpG  methylation  at  the  promoter  of  PHD3  negative 
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Figure  3.  DNA  methylation  analysis  of  the  PHD3  5'-CpG  island  in  human  melanoma,  prostate  and  breast  carcinoma  cell  lines.  A) 

Schematic  representation  of  the  CpG  distribution  in  the  5'-regulatory  region  of  the  PHD3  gene  from  UCSC  Genome  Browser  coordinates 
bp3442051 1-bp34419838.  The  CpG  sites  are  represented  by  vertical  tick  marks,  and  the  beginning  of  exon  1  is  depicted  as  a  bent  arrow.  The 
locations  of  PCR  primers  used  to  assess  genomic  methylation  (SB3  and  SB4)  and  chromatin  accessibility  (CA1  and  CA2),  respectively.  B)  Bisulfite 
sequencing  was  performed  on  DNA  isolated  from  two  non-tumorigenic  epithelial  cell  lines  (NPrEC  and  hTERT-HMEI)  and  six  tumor  cell  lines  (DU  145, 
22RV.1,  PC-3,  MB-435,  HS578T,  and  MCF7).  Each  grid  contains  at  least  5  rows  representing  the  5  individually  cloned  and  sequenced  bisulfite  PCR 
products  from  the  indicated  cell  lines.  Each  row  contains  58  boxes  representing  the  58  CpGs  in  an  individual  cloned  molecule.  Open  and  filled  boxes 
are  unmethylated  and  methylated  CpG  sites,  respectively. 
doi:1 0.1 371  /journal.pone.001 461 7.g003 


ACHN  cells,  suggesting  an  alternative  mechanism  for  silencing  in 
this  cell  line. 

PHD3  promoter  methylation  status  does  not  correlate 
with  hypoxia  induced  HIF-la  protein  accumulation  or  HIF 
transcriptional  activity 

The  presence  of  PHD3  promoter  methylation  in  such  a  broad 
range  of  epithelial  malignancies  suggests  that  it  may  be  a  selective 
advantage  for  tumor  survival.  One  hypothesis  is  that  PHD3  silencing 
by  promoter  methylation  may  allow  for  an  increased  HIF 
transcriptional  response  during  hypoxic  conditions.  In  order  to 
determine  whether  PHD3  promoter  methylation  specifically  affects 
the  hypoxia  response  pathway,  we  performed  western  blots  on  cell 
lysates  front  MCF7,  PC-3,  MB-435  and  DU  145  cell  lines  to  compare 
the  HIF  protein  levels  and  the  HIF  transcriptional  response  to 
hypoxia  (Figure  5a).  Following  24  hours  of  hypoxia,  HIF-la 
protein  was  upregulated  in  all  the  cell  lines  regardless  of  PHD3 
expression  status.  We  also  observed  that  DU  145  cells  appear  not  to 


express  detectable  levels  of  PHD3  protein.  We  are  unsure  whether 
this  is  due  to  limits  of  detection  by  our  PHD  3  antibody,  or  if  DU  145 
cells  downregulate  PHD3  expression  by  a  posttranslational  mecha¬ 
nism.  Our  antibody  appears  to  be  specific  to  PHD3  as  transduction  of 
an  adenoviral-PHD3  expression  vector  into  PC3  cells  produces  a 
band  at  an  identical  molecular  weight  as  the  band  seen  in  MCF7  cells 
(Figure  S3).  Interestingly  though,  MCF7  cells,  which  do  express 
basal  levels  of  PHD  3  mRNA  and  protein,  displayed  the  largest 
induction  of  HIF- 1  a  protein.  Thus,  HIF- 1  a  protein  levels,  in  general, 
did  not  show  any  correlation  with  presence  or  absence  of  PHDS.  We 
also  found  HIF-2a  to  be  expressed  under  normoxic  conditions  in 
MCF7.  Moreover,  MB-435  cells,  which  express  the  lowest  levels  of 
PHD3  niRNA  out  of  all  the  cell  lines  tested  appear  not  to  express 
HIF-2a  at  an  appreciable  level.  Thus,  loss  of  PHD 3  does  not  appear 
to  be  significantly  correlated  with  an  accumulation  ofHIF-la  or  HIF- 
2a  levels  in  these  cell  lines. 

To  further  investigate  the  effect  of  PHDS  promoter  methylation 
on  the  transcriptional  response  of  the  hypoxia  response  pathway, 
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Figure  4.  The  methylated  PHD3  gene  in  melanoma,  prostate,  breast  and  renal  carcinoma  cell  lines  is  refractory  to  induction  by 
hypoxia.  A)  Melanoma,  prostate  and  breast  carcinoma  cell  lines  were  treated  with  hypoxia  (1%  02)  or  normoxia  (21%  02)  for  24  hours.  Total  RNA 
was  extracted  and  converted  to  cDNA  by  reverse  transcription.  Quantitative  real-time  reverse  transcription-PCR  analysis  of  PHD3  was  performed  with 
normalization  to  GAPDH  gene  expression.  Relative  quantitation  was  determined  by  the  DDCt  method.  ND  =  not  detectable.  Error  bars  =  SEM.  n  =  3. 
B)  Renal  clear  cell  carcinoma  cell  lines  were  either  untreated  or  treated  with  hypoxia  as  in  A).  The  right  panel  depicts  relative  PHD3  mRNA  levels 
compared  to  mRNA  extracted  from  normal  renal  tissue.  The  left  panel  depicts  the  methylation  status  at  each  of  58  CpG  dinucleotides  present  in  the 
PHD3  CpG  island  of  the  representative  renal  carcinoma  cell  lines.  Error  bars  =  SEM.  n  =  3. 
doi:1 0.1 371  /journal.pone.001 461 7.g004 


we  transfected  MCF7,  PC-3,  MB-435  and  DU  145  cancer  cell 
lines  with  an  HRE-luciferase  reporter  construct  [19].  Following 
24  hours  of  hypoxia,  luciferase  activity  was  measured  and  plotted 
relative  to  luciferase  activity  in  normoxic  cells  (Figure  5B).  There 
was  no  correlation  between  PHD3  promoter  methylation  status 
and  hypoxic  induction  of  luciferase.  The  results  of  this  experiment 
resembled  the  pattern  of  HIF-la  protein  accumulation  seen  in 
figure  5A.  MCF7  cells  showed  the  largest  induction  luciferase 
activity,  whereas  the  other  cell  lines,  which  do  not  express 
detectable  levels  of  PHD3  protein,  were  comparable  to  one 
another.  Although  not  absolute,  these  data  are  highly  suggestive 
that  F1IF  protein  stabilization  and  transcriptional  activity  is  largely 
independent  of  PHD3  expression. 


PHD3  promoter  methylation  is  absent  in  primary  human 
prostate  adenocarcinomas 

Recent  data  published  by  Flatzimichael  et  al.  suggested  that 
acquisition  of  PHD3  promoter  methylation  may  be  a  relatively 
common  event  in  certain  plasma  cell  neoplasias  [17].  There¬ 
fore,  we  asked  whether  primary  human  prostate  neoplasias 
contained  methylation  at  the  PHD3  CpG  island.  We  extracted 
DNA  and  RNA  from  frozen  sections  of  10  prostate  cancer 
specimens  containing  a  minimum  of  70%  malignant  tissue  with 
a  Gleason  score  ranging  from  7-9  as  well  as  3  benign  prostate 
specimens.  Real-time  PCR  of  extracted  RNA  showed  that  all 
10  tumors  contained  decreased  PFID3  mRNA  expression 
compared  to  3  benign  tissue  specimens  (Figure  6A  and  B). 
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Figure  5.  The  HIF-1  response  is  not  impaired  in  prostate  and  breast  carcinoma  cell  lines  that  lack  PHD3.  A)  Melanoma,  breast  and 
prostate  cancer  cell  lines  with  methylated  PHD3  promoters  (MB-435,  PC-3)  and  non-methylated  promoters  (MCF7,  DU  145)  were  subjected  to  hypoxia 
(1%  02)  or  normoxia  for  24  hours.  Thirty  micrograms  of  whole  cell  lysate  was  western  blotted  for  the  presence  of  HIF-1  or,  HIF-2oc  and  PHD3;  actin  was 
used  as  a  loading  control.  B)  Breast  and  prostate  cancer  cell  lines  were  transfected  with  an  HRE-luciferase  reporter  vector  and  Renilla  luciferase  vector 
and  then  subjected  to  hypoxia  (1%  02)  or  normoxia  for  24  hours.  Luciferase  activities  after  24  hours  of  hypoxia  were  determined  and  are  depicted 
relative  to  luciferase  activities  in  cells  under  normoxia.  Error  bars  =  SEM.  n  =  3  for  MCF7  and  PC-3.  n  =  2  for  MB-435  and  DU  145. 
doi:1 0.1 371  /journal. pone.001 461 7.g005 


Probing  for  methylated  CpGs  in  select  tumor  specimens  using 
the  MethylMiner  kit  indicated  the  possible  presence  of 
methylated  CpGs  at  the  PHD3  CpG  island  in  tumor  sample  1 
(Figure  S4),  which  also  contained  the  lowest  PHD3  mRNA 
levels.  However,  further  bisulfite  sequencing  of  8  clones  from 


tumor  sample  1  and  tumor  sample  4  did  not  detect  any 
methylated  CpGs  (data  not  shown).  This  could  be  due  to  better 
sensitivity  to  population  average  methylation  states  as  com¬ 
pared  to  the  single  molecule  at  a  time  approach  of  bisulfite 
sequencing. 
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Figure  6.  PHD3  mRNA  expression  is  downregulated  in  multiple  primary  human  prostate  cancer  specimens.  A)  Total  mRNA  was 
isolated  from  frozen  sections  of  primary  human  prostate  cancer  specimens  with  Gleason  scores  ranging  from  7-9.  Quantitative  real-time  PCR  was 
performed  using  PHD3  specific  TaqMan  primer-probe.  Relative  PHD3  mRNA  expression  for  each  tumor  sample  is  represented  as  (average  dCt  of  n  =  3 
benign  prostate  tissue  samples)  -  (dCt  tumor  sample).  Samples  were  normalized  to  GAPDH.  Dotted  lines  represent  +/—  1  SD  for  benign  tissue  PHD3 
mRNA  expression.  B)  Box  plot  depicting  PHD3  mRNA  expression  from  samples  shown  in  A.  p  value  is  based  on  ANOVA  between  3  benign  samples 
and  10  malignant  samples. 
doi:1 0.1 371  /journal.pone.001 461 7.g006 


Discussion 

Perturbations  in  the  cellular  responses  to  hypoxia  are  well  known 
to  play  a  role  in  the  malignant  process.  Familial  mutations  in  VHL,  a 
negative  regulator  of  the  HIF-a  proteins,  results  in  vascular  tumors 
of  the  brain,  spinal  cord  and  retina,  as  well  as  appearance  of  renal 
clear-cell  carcinomas  [23].  PHD  proteins  play  a  role  upstream  of 
VHL  regulation;  they  hydroxylate  HIF-a  proteins,  creating  a 
binding  site  for  VHL  [24,25,26].  Thus,  it  is  feasible  that 


deregulation  of  PHD  activity  or  expression  could  also  contribute 
to  the  malignant  process.  In  fact,  an  absence  of  PHD3  upregulation 
following  hypoxia  has  been  observed  in  multiple  human  cell  lines 
from  tumors  of  the  breast,  prostate  and  brain  [10,16].  A  recent 
clinical  study  of  breast  tumors  containing  BRCA  mutations  supports 
the  hypothesis  that  PHD3  plays  an  important  role  in  malignancy. 
This  study  found  a  positive  correlation  between  decreased  PHD3 
expression  and  a  basal  phenotype,  which  is  considered  a  higher 
grade  and  more  aggressive  tumor  [27]. 
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Here,  we  report  aberrantly  silenced  basal  mRNA  expression  of 
PHD3  in  breast,  prostate,  melanoma  and  renal  cell  carcinoma  cell 
lines,  and  the  absence  of  PHD3  mRNA  induction  upon  hypoxic 
stimulus.  PHD3  expression  could  be  recapitulated  in  some  PHD3 
negative  cell  lines  after  treatment  with  5-Aza-dC,  a  DNA 
methyltransferase  inhibitor,  implicating  DNA  methylation  as  a 
mechanism  for  the  decreased  expression  of  PHD  3  mRNA  in  these 
cell  lines.  PHD3  promoter  methylation  was  verified  by  sequence 
analysis  of  PGR  products  cloned  from  bisulfite-treated  genomic 
DNA.  We  found  that  among  the  human  cancer  cell  lines 
investigated,  PC-3,  MB-435,  HS578T,  and  769-P  cell  lines  have 
hypermethylated  PHD3  CpG  islands.  Furthermore,  the  PHD3 
promoter  region  was  more  resistant  to  DNase  I  in  PC-3  cells 
(hypermethylated  PHD3  promoter)  compared  to  MCF7  cells 
(hypomethylated  PHD3  promoter).  The  methylation  of  the  PHD3 
promoter  in  these  carcinoma  cell  lines  appears  to  be  aberrant  since 
insignificant  DNA  methylation  was  found  in  the  non-transformed 
cell  counterparts  of  prostate  and  mammary  epithelial  cell  lines 
NPrEC  and  HME1  respectively.  The  apparently  aberrant  PHD3 
promoter  methylation  status  in  these  cell  lines  is  the  likely 
mechanism,  at  least  in  part,  for  PHD3  transcriptional  suppression 
because  DNA  methylation  is  typically  associated  with  a  condensed 
heterochromatin  state,  and  is  known  to  inhibit  transcription  factor 
binding  to  promoter  regions  of  genes  [28] . 

Although  all  members  of  the  PHD  family  have  the  ability  to 
hydroxylate  both  HIF-1  and  HIF-2a,  the  specificities  appear  to 
differ  slightly.  PHD2  has  been  reported  to  play  a  more 
pronounced  role  in  the  regulation  of  HIF-1  a,  whereas  PHD  3 
more  strongly  affects  HIF-2a  stability  [10].  Therefore,  it  would 
seem  likely  that  loss  of  PHD3  expression  by  promoter  methylation 
would  convey  a  cellular  advantage  mediated  through  increased 
HIF-la  and/or  HIF-2a  stability  during  hypoxia.  This  could  lead 
to  increased  expression  of  VEGF  and  erythropoietin  with 
subsequent  vascular  recruitment.  Tumors  of  the  breast,  skin, 
kidney  and  prostate,  being  solid  tumors,  would  certainly  benefit 
from  an  increase  in  vascular  supply  to  hypoxic  areas.  In  fact,  diere 
is  evidence  that  cell  lines  from  other  solid  tumors  downregulate 
PHD3  as  well.  Henze  et  al.  have  shown  that  several  glioma  cell 
lines  display  little  to  no  PHD  3  protein  expression  during 
normoxia,  with  no  induction  upon  hypoxia  when  compared  to 
other  glioma  cell  lines  studied  [16].  Our  results  would  predict  that 
a  subset  of  those  cell  lines  have  aberrant  methylation  of  the  PHD3 
CpG  island.  Interestingly,  HIF-2ct  protein  levels  after  18  hours  of 
hypoxia  appeared  lower  in  PHD3  non-expressing  cells  than  those 
in  cells  expressing  relatively  high  levels  of  PHD3  [16].  This  is  the 
opposite  of  what  we  had  expected,  and  may  demonstrate  the 
ability  of  other  PHD  family  members  to  substitute  for  the  loss  of 
PHD  3  expression  in  regulating  PIIF-la  and  HIF-2a  stability.  This 
hypothesis  is  supported  by  data  front  Appelhoff  et  al.,  who 
measured  relative  protein  amounts  of  PHD  isoforms  in  multiple 
cell  lines.  In  cell  MB-435  and  ZR751  cell  lines  where  PHD3  is  low 
or  absent,  PFID2  protein  is  greatly  elevated.  [29] .  Similarly,  we 
found  that  MCF7  cells,  which  basally  express  PHD3  at  the  ntRNA 
and  protein  level  during  normoxic  conditions,  also  express 
relatively  high  levels  of  HIF-2a.  Taken  together,  we  interpret 
these  findings  to  suggest  that  PHD3  silencing  by  CpG  methylation 
may  not  have  a  significant  impact  on  HIF-1  and  HIF-2a  protein 
levels  in  the  cell  lines  that  we  tested. 

In  fact,  the  results  of  our  study  support  a  mechanism  whereby 
PHD  3  silencing  by  PHD3  CpG  island  methylation  affects 
pathways  outside  of  the  conventional  hypoxic  response  pathway. 
When  HIF  transcriptional  activity  was  measured  through  a 
hypoxia  responsive  HRE-luciferase  reporter,  we  observed  nearly 
equal  transcriptional  responses  to  hypoxia  in  3  out  of  4  cell  lines, 


which  included  both  PHD3  silenced  as  well  as  PHD3  expressing 
cell  lines.  In  MCF7  cells,  which  express  PHD3  at  the  protein  level 
and  do  not  contain  PHD3  promoter  methylation,  we  observed 
nearly  a  15-fold  increase  in  luciferase  induction  upon  hypoxic 
treatment.  These  results  suggest  that  PHD3  expression  status  does 
not  significantly  affect  HIF  protein  stabilization  or  HIF  transcrip¬ 
tional  activity  through  an  HRE  containing  promoter  upon 
exposure  to  1  %  oxygen.  Therefore,  modulation  of  an  alternative 
cellular  pathway(s)  remains  an  open  candidate  for  mediating  the 
effects  of  PHD3  loss  in  malignancies. 

Besides  the  HIF  family  proteins,  other  interacting  partners  of 
PHD3  have  already  been  discovered.  PHD3  is  a  known  player  in 
both  neuronal  apoptosis  and  in  myoblast  differentiation  [30,31] 
PHD3  also  appears  to  interact  with  Bcl-2  to  induce  apoptosis  in 
H9c2  cells  in  response  to  doxyrubicin  [32]. Furthermore,  PHD3 
has  also  been  reported  to  destabilize  ATF-4  through  a  novel 
oxygen-dependant  domain  on  ATF-4  [33].  ATF-4  is  involved  in 
the  regulation  of  angiogenesis  and  metabolism  [34],  Thus, 
upregulation  of  ATF-4  by  a  loss  of  PHD3  could  promote  cell 
survival.  In  addition,  PHD  3  has  also  been  recently  reported  to 
inhibit  IKK(3.  An  increase  in  IKKfl  activity  in  the  absence  of 
PHD  3  could  confer  a  growth  advantage  to  cells  through  an 
increase  NFkB  signaling  [35].  Constitutive  NFkB  activity  is  an 
important  and  common  event  in  T-  and  B-cell  derived 
malignancies  [36,37]  and  may  explain  the  PHD3  promoter 
methylation  recently  reported  in  plasma  and  B-cell  neoplasia  by 
Hatzimichael  et  al.  [17].  These  interesting  prospects  related  to 
identification  of  the  downstream  targets  of  PHD  3  signaling  will 
undoubtedly  become  a  focus  of  future  investigations. 

Our  inability  to  detect  PHD3  promoter  DNA  methylation  in 
primary  human  prostate  tumors  was  surprising,  however  our 
results  are  supported  by  a  recent  study  by  Huang  et  al.  [38],  who 
screened  168  invasive  breast  carcinomas  and  did  not  find  evidence 
of  PPID3  DNA  methylation  using  melting  curve  analysis  of  bisulfite 
converted  DNA.  We  cannot  rule  out  the  possibility  that  a  low  level 
of  semi-methylated  PHD3  CpG  islands  are  present  in  their  samples 
as  the  “intermediately  methylated”  controls  were  not  truly 
intermediately  methylated,  but  rather  were  mixtures  of  100% 
methylated  DNA  with  0%  methylated  DNA.  Nonetheless,  it  seems 
clear  that  a  large  proportion  of  primary  epithelial  tumors  does  not 
contain  a  high  degree  of  PHD3  promoter  methylation,  and  may 
not  be  the  ideal  specimens  for  detection  of  methylation  at  this 
locus.  Given  PHD3's  purported  ability  to  negatively  regulate  the 
NFkB  pathway  [30,35],  and  the  widely  reported  involvement  of 
NFkB  in  cell  migration  and  metastasis  [39,40],  future  studies  on 
PHD3  CpG  island  methylation  in  clinical  samples  of  metastatic 
disease  as  opposed  to  primary  tumors  may  yield  more  positive 
results. 

Here,  we  are  the  first  to  report  DNA  methylation  of  the  PHD3 
CpG  island  in  solid  tumor  cell  lines  derived  from  diverse  cell  types. 
PHD3  methylation  in  carcinoma  cells  was  associated  with  their 
inability  to  appropriately  upregulate  PHD3  mRNA  upon  exposure 
to  hypoxia.  We  are  also  the  first  to  show  evidence  that  this 
aberrant  expression  fails  to  correlate  with  an  increase  in  HIF 
protein  accumulation  and  transcriptional  activity  upon  exposure 
to  hypoxia  in  the  cell  lines  examined.  The  presence  of  PHD3 
promoter  hypermethylation  and  PHD  3  silencing  in  such  a  wide 
range  of  cancer  types  suggests  this  might  be  a  common  event  that 
elicits  a  selective  advantage  for  tumors.  Our  data  suggest  that  at 
least  in  some  cell  lines,  the  nature  of  this  advantage  may  extend 
beyond  hypoxia  resistance.  Furthermore,  the  selective  event  may 
occur  during  or  after  the  process  of  invasion/ metastasis,  as  we  and 
others  have  not  found  evidence  of  methylation  in  primaiy  solid 
tumors  [38], 
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Supporting  Information 

Figure  SI  Validation  of  methylated  DNA  enrichment  as  a  tool 
for  detecting  methylated  CpG  regions  in  human  genomic  DNA. 
A)  Genomic  DNA  from  human  melanoma  and  breast  cancer  lines 
was  enriched  for  methylated  CpG  dinucleotides  at  the  PHD3  CpG 
island  using  the  MethylMiner  kit,  followed  by  quantitative  real¬ 
time  PCR  analysis  using  PHD3  CpG  island-specific  PCR  primers. 
S  =  supernatant,  representing  unmethylated  DNA.  El  -  E5 
represents  elutions  with  200,  650,  1100,  1550,  2000  mM  NaCl 
respectively.  Input  =  l/60th  total  input  DNA.  Total  amounts  of 
eluted  DNA  from  each  fraction  are  represented  as  a  fraction  of 
input  (left).  Methylated  CpG  sequences  from  bisulfite-converted 
DNA  in  the  corresponding  cell  lines  are  depicted  for  comparison 
(right).  B)  Control  100%  methylated  and  0%  methylated  oligos 
supplied  by  the  MethylMiner  kit  were  subjected  to  the 
MethylMiner  protocol.  S  =  supernatant,  representing  unmethy¬ 
lated  DNA.  E1-E5  represent  elutions  with  200,  650,  1100,  1550, 
2000  mM  NaCl  respectively.  Identical  volumes  of  DNA  from  each 
eluate  were  subjected  to  quantitative  real-time  PCR  using  primers 
supplied  by  the  MethylMiner  kit.  DNA  content  in  each  fraction  is 
represented  as  arbitrary  units. 

Found  at:  doi:10.1371/journal.pone. 0014617. sOOl  (10.36  MB 

TIF) 

Figure  S2  The  methylated  PHD3  gene  in  non-expressing  cells  is 
maintained  in  a  less  accessible  state  than  the  non-nrethylated 
PHD3  gene  in  expressing  cells.  A)  Nuclei  from  PHD3-positive 
MCF7  and  PHD3-negative  PC-3  carcinoma  cell  lines  were 
isolated  and  enzymatically  restricted  with  DNase  I.  Primers  CA1 
and  CA2  (see  Fig.  3A)  were  used  for  quantitative  real-time  PCR 
(right  panels)  to  amplify  a  region  also  assessed  for  cytosine 
methylation.  Accessibility  indices  (left  panels)  were  calculated  as 
follows:  AI  =  2((Ct  DNase  treated)  -  (Ct  Untreated)).  B)  GAPDH 
accessibility  indices  were  simultaneously  assessed  as  a  control  for  a 
constitutively  expressed  gene  in  both  cell  lines. 

Found  at:  doi:10.1371/journal.pone.0014617.s002  (3.00  MB  TIF) 

Figure  S3  PHD3  antibody  specificity.  PC3  cells  were  transduced 
with  an  increasing  MOI  of  adenoviral-PHD3  vector.  Western  blot 
using  Novusl00-139  antibody  co-incubated  with  anti-actin 
antibody  indicated  an  approximately  27  kDa  band  in  MCF7  cells 
that  migrates  at  the  same  molecular  weight  as  a  band  present  in 
PHD3  transduced  PC3  cells. 

Found  at:  doi:10.1371/journal.pone. 0014617. s003  (3.00  MB 
TIF) 
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Reactive  oxygen  species  (ROS)  are  critical  in  a  broad  spectrum  of  cellular  processes  including  signaling,  tumor 
progression,  and  innate  immunity.  The  essential  nature  of  ROS  signaling  in  the  immune  systems  of  Drosophila 
and  zebrafish  has  been  demonstrated;  however,  the  role  of  ROS,  if  any,  in  mammalian  adaptive  immune 
system  development  and  function  remains  unknown.  This  work  provides  the  first  clear  demonstration  that 
thymus-specific  elevation  of  mitochondrial  superoxide  (02_)  disrupts  normal  T  cell  development  and  impairs 
the  function  of  the  mammalian  adaptive  immune  system.  To  assess  the  effect  of  elevated  mitochondrial 
superoxide  in  the  developing  thymus,  we  used  a  T-cell-specific  knockout  of  manganese  superoxide  dismutase 
(i.e.,  SOD2)  and  have  thus  established  a  murine  model  to  examine  the  role  of  mitochondrial  superoxide  in  T 
cell  development.  Conditional  loss  of  SOD2  led  to  increased  superoxide,  apoptosis,  and  developmental  defects 
in  the  T  cell  population,  resulting  in  immunodeficiency  and  susceptibility  to  the  influenza  A  virus  H1N1.  This 
phenotype  was  rescued  with  mitochondrially  targeted  superoxide-scavenging  drugs.  These  findings 
demonstrate  that  loss  of  regulated  levels  of  mitochondrial  superoxide  lead  to  aberrant  T  cell  development 
and  function,  and  further  suggest  that  manipulations  of  mitochondrial  superoxide  levels  may  significantly 
alter  clinical  outcomes  resulting  from  viral  infection. 

©  2010  Elsevier  Inc.  All  rights  reserved. 


The  term  reactive  oxygen  species  (ROS)  describes  a  large  group  of 
free  radical  and  non-free  radical  oxygen-containing  compounds  (e.g., 
superoxide,  02~ ;  hydrogen  peroxide,  H202;  peroxynitrite,  ONOO~; 
hydroxyl  radical,  'OH)  [1],  It  is  commonly  accepted  that  ROS  are  by¬ 
products  of  normal  metabolism  and  as  such  act  to  damage  cellular 
components  such  as  nucleic  acids,  proteins,  or  lipids  [2-4].  Because  of 
this,  ROS  have  been  implicated  in  many  different  diseases  including 
cancer,  atherosclerosis,  amyotrophic  lateral  sclerosis,  Alzheimer 
disease,  and  many  others  [5-8].  One  specific  role  of  ROS  is  their 
ability  to  enhance  the  pathogenesis  of  infections,  such  as  influenza 


Abbreviations:  ROS,  reactive  oxygen  species;  0/,  ,  superoxide;  SOD2,  manganese 
superoxide  dismutase;  1AV,  influenza  A  virus;  ONOO  ,  peroxynitrite;  "OH,  hydroxyl 
radical;  SOD1,  copper/zinc  superoxide  dismutase;  SOD3,  extracellular  superoxide 
dismutase;  SOD2  /  ,  conditional  SOD2  knockout  mouse;  S0D2L/L,  control  loxP- 
containing  mouse;  DHE,  dihydroethidium;  DCFH-DA,  dihydrodichlorofluorescein 
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[9,10],  It  has  been  demonstrated  that  during  times  of  influenza 
infection  ROS  may  damage  lung  parenchyma  cells,  but  that  this  injury 
may  be  ameliorated  by  antioxidant  supplementation  [11-14].  Current 
theories  propose  the  mechanism  behind  this  benefit  to  be  attenuation 
of  ROS  produced  by  the  innate  immune  system,  but  this  is  not 
commonly  accepted  and  is  still  highly  debated.  More  recently  it  has 
been  shown  that  cells  possess  the  ability  to  exploit  ROS  for  signaling 
and  functional  purposes.  For  example,  many  transcription  factor 
pathways  are  sensitive  to  oxidative  stress  and  as  such  are  able  to  help 
cells  adapt  to  large  deviations  in  redox  status  [15-18].  Moreover,  ROS 
are  essential  in  the  development  of  certain  organ  systems  and  even 
whole  organisms  [19,20].  With  this  knowledge,  the  importance  of  ROS 
in  biology  is  being  elucidated,  but  many  questions  about  tissue- 
specific  dependence,  specific  ROS  functions,  molecular  targets  of  ROS, 
and  ROS  mechanisms  of  action  remain  unanswered. 

One  organ  system  in  which  ROS  have  been  widely  described  is 
the  immune  system.  The  biological  relevance  of  ROS  in  this  system 
was  first  depicted  when  it  was  found  that  leukocytes  generated 
ROS  during  the  oxidative  burst  to  neutralize  pathogens  [21].  More 
recent  studies  have  demonstrated  the  importance  of  ROS  in  the 
downstream  intracellular  signaling  post-T-cell  activation  [22-24]. 
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Furthermore,  evidence  has  shown  that  hydrogen  peroxide  acts  as 
an  important  chemoattractant  to  direct  leukocytes  to  wound  margins 
at  sites  of  injury,  which  was  pivotal  in  demonstrating  that  immune 
cells  are  able  to  respond  to  exogenous  ROS  in  addition  to  producing 
endogenous  levels  [25],  Recent  studies  have  demonstrated  the 
role  of  ROS  in  priming  the  development  of  the  primitive  immune 
system  in  Drosophila,  illuminating  the  first  described  function  of  ROS 
in  the  development  of  this  organ  system  in  an  invertebrate  species 
[26],  Finally,  it  has  been  demonstrated  that  an  intracellular  pro¬ 
oxidant  shift  occurs  before  normal  thymocyte  apoptosis  during 
selection  [27],  and  numerous  groups  have  shown  increased  ROS 
levels  in  pathologic  T  cell  systems  [28-30],  In  contrast,  few  studies 
to  date  have  focused  on  how  pro-  or  antioxidants  affect  the  normal 
mammalian  adaptive  immune  system,  and  to  our  knowledge  no 
published  studies  have  directly  addressed  the  role  of  excess 
mitochondrial  superoxide  in  the  growth  or  function  of  the  adaptive 
immune  system.  With  the  understanding  that  ROS  play  a  major  part 
in  intracellular  signaling  and  cellular  damage,  we  hypothesized  that 
mitochondrial  superoxide  may  play  a  distinct  role  in  normal 
development  and  maintenance  of  the  mammalian  adaptive  immune 
system,  which  could  be  central  to  intercellular  communication 
between  the  innate  and  the  adaptive  branches  of  the  immune 
system. 

Although  the  conceptual  framework  underlying  this  hypothesis 
seems  simple  in  nature,  the  availability  of  animal  models  with  which 
to  address  this  hypothesis  is  limited.  Numerous  constitutive  antiox¬ 
idant  enzyme  knockout  animals  have  provided  valuable  information 
on  the  global  developmental  importance  of  these  proteins  [31-36], 
but  these  studies  have  not  assessed  the  role  of  specific  ROS  in  tissue- 
specific  adaptive  immunity.  To  address  these  deficits,  we  used  a 
conditional  T  cell  manganese  superoxide  dismutase  (i.e.,  SOD2) 
knockout  mouse  to  examine  the  role  of  increased  steady-state  levels 
of  superoxide  during  mammalian  adaptive  immune  system  develop¬ 
ment  [37].  The  superoxide  dismutase  class  of  enzymes  specifically 
scavenges  superoxide  in  biological  systems  [38].  Mammals  contain 
three  variants  of  the  enzymes:  cytoplasmic  Cu/Zn  SOD,  SOD1; 
mitochondrial  SOD2;  and  extracellular  SOD,  SOD3.  Because  SOD2 
has  the  explicit  role  of  eliminating  mitochondrial  superoxide,  tissue- 
specific  disruption  of  this  activity  should  provide  an  excellent  in  vivo 
model  under  conditions  where  compartmentalized  superoxide 
metabolism  is  disrupted.  A  constitutive  SOD2  knockout  mouse  has 
been  created,  but  because  of  the  postnatal  developmental  dependence 
on  SOD2  the  animal  succumbs  to  numerous  organ  failures  shortly 
after  birth  [39,40].  Furthermore,  because  of  the  mouse's  limited  life 
span  no  examination  of  the  immune  system  was  reported.  Taken 
together,  our  model  serves  as  the  first  described  animal  model  for 
studying  the  effects  of  perturbing  steady-state  mitochondrial  super¬ 
oxide  levels  on  the  development  and  function  of  the  mammalian  T  cell 
adaptive  immune  system. 

Materials  and  methods 

Mice 

Mice  homozygous  for  the  floxed  SOD2  allele  (i.e.,  B6.Cg -Sod2tml, 
shorthand  SOD2L/L),  in  which  exon  3  of  the  SOD2  gene  is  flanked  by 
two  loxP  sequences,  have  been  previously  described  [37],  B6.Cg-Tg- 
Lck-Cre548Jxm/},  or  Lck-Cre,  mice  (in  which  Cre  recombinase  is 
exogenously  expressed  under  the  control  of  the  proximal  lympho¬ 
cyte-specific  kinase  and  Cre  recombinase  becomes  activated  in  a(3  T 
cells  during  the  CD4~/CD8~  to  CD4+/CD8+  stage  of  development) 
and  B6.Cg-Tg-\7avJ -iCre/>2K'7J,  or  Vav-iCre,  mice  (Cre  recombinase  is 
exogenously  expressed  under  control  of  the  Vav  promoter,  and  Cre 
recombinase  becomes  activated  within  the  hematopoietic  stem  cell  and 
affects  all  lymphoid,  myeloid,  and  erythroid  lineages)  were  generously 
donated  by  Dr.  Adam  Dupuy  (The  University  of  Iowa)  and  have  been 


previously  described  [41,42].  B6.Cg-Tg -AlbCrez,Mgn /]  mice  (Cre  recom¬ 
binase  is  exogenously  expressed  under  the  control  of  the  albumin 
promoter  and  as  such  expression  is  limited  to  the  liver)  were  generously 
donated  by  Dr.  Curt  Sigmund  (The  University  of  Iowa)  and  have  been 
previously  described  [43].  Finally,  B6.129-Tg-MM7VCre4Mam/J  mice  (Cre 
recombinase  is  exogenously  expressed  under  the  control  of  the  mouse 
mammary  tumor  virus  promoter,  and  expression  is  primarily  limited  to 
mammary  tissue)  were  purchased  from  The  Jackson  Laboratory  and 
have  been  previously  described  [44].To  obtain  conditional  T  cell  SOD2 
homozygous  knockout  animals  (i.e.,  SOD2_/_),  parent  strains  of  both 
the  floxed  SOD2  and  the  Lck-Cre  mice  were  bred  to  generate  FI 
heterozygotes  (i.e.,  SOD2wt/_).  The  FI  generation  was  then  bred  back  to 
the  parent  floxed  SOD2  mice  to  create  F2  homozygous  knockouts.  Lck- 
Cre  was  passed  only  through  male  parents  to  limit  nonspecific  oocyte 
expression.  The  mice  used  were  of  pure  C57BL/6  background,  and 
littermate  floxed  animals  ( i.e.,  SOD2L/L)  served  as  controls.  All  work  was 
performed  under  the  approval  of  the  Institutional  Animal  Care  and  Use 
Committee  at  The  University  of  Iowa. 

Tissue  isolation 

For  all  experiments,  fresh  tissue  was  harvested  and  used  from 
animals  at  6  weeks  of  age  unless  otherwise  noted.  For  thymocyte, 
splenocyte,  and  lymph  node  preps,  organs  were  freshly  harvested 
upon  necropsy  and  placed  in  a  solution  of  Hanks’  buffered  salt 
solution  (HBSS)  containing  10%  fetal  bovine  serum.  Tissues  were 
dissociated  by  physical  disruption  using  ground  glass.  After  this, 
lymphocytes  were  isolated  by  mouse  optimized  FicoLite-LM  gradient 
(Atlanta  Biologicals,  Atlanta,  GA,  USA).  Remaining  red  blood  cell 
contaminants  were  removed  by  red  blood  cell  lysis  buffer,  and  the 
remaining  lymphocytes  were  washed  twice  with  PBS  before  further 
analysis.  Peripheral  blood  was  isolated  by  retro-orbital  bleeding  into 
heparin-coated  capillary  tubes. 

Real-time  PCR 

RNA  was  extracted  using  the  Trizol  method  and  was  quantified  by 
the  use  of  a  Nanodrop  ND-1000.  One  microgram  from  each  sample 
was  reverse  transcribed  using  the  ABI  cDNA  archive  kit.  Generated 
cDNA  was  then  subjected  to  SYBR  green  quantitative  real-time  PCR 
with  primers  specific  to  each  individual  transcript  and  the  18  S  control 
(supplementary  primer  sequences)  under  the  following  PCR  para¬ 
meters:  95  °C  for  10  min,  followed  by  35  cycles  of  95  °C  for  15  s  and 
60  °C  for  1  min.  A  threshold  in  the  linear  range  of  PCR  amplification 
was  selected  and  the  cycle  threshold  (Ct)  determined.  Levels  of 
transcripts  were  then  normalized  to  the  18  S  control  and  compared 
relative  to  the  SOD2L/L  control  using  the  AACt  method. 

Western  blot  analysis 

Protein  was  extracted  using  standard  RIPA  buffer  and  was 
quantified  on  the  basis  of  the  Bradford  assay  and  a  standard  curve. 
Protein  was  run  on  an  SDS-PAGE  gel  for  separation  of  different-sized 
proteins.  The  products  were  then  transferred  from  the  gel  to  a 
nitrocellulose  membrane.  Identification  and  quantification  of  the 
amount  of  protein  were  performed  by  addition  of  specific  primary 
antibody  and  then  the  addition  of  a  horseradish  peroxidase-tagged 
secondary  antibody  to  the  first.  Chemiluminescent  substrate  was 
added  to  the  blot  for  exposure  of  products,  and  film  was  exposed  and 
qualitatively  quantified  for  protein  amount.  Loading  errors  were 
controlled  for  by  normalizing  to  [3-actin.  Antibodies  used  for  this 
study  were  SOD2  (Millipore,  Billerica,  MA,  USA),  aconitase  (Abeam, 
Cambridge,  MA,  USA),  SDHB  (Abeam),  LC3  (Cell  Signaling  Technology, 
Beverly,  MA,  USA),  p53  (Abeam),  and  p-actin  (Abeam). 
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SOD2  activity 

To  analyze  the  activity  of  the  superoxide  dismutase  enzymes,  the 
indirect  competitive  inhibition  assay  developed  by  Spitz  and  Oberley 
was  used  [45].  Briefly,  superoxide  is  generated  from  xanthine  by 
xanthine  oxidase  and  detected  by  recording  the  reduction  of  nitroblue 
tetrazolium  (NBT).  SOD  scavenges  superoxide  and  competitively 
inhibits  the  reduction  of  NBT.  One  unit  of  SOD  activity  is  defined  as  the 
amount  of  protein  required  to  inhibit  50%  of  the  maximal  NBT 
reduction.  To  obtain  the  amount  of  MnSOD  activity,  sodium  cyanide 
(5  mmol/L)  is  added  to  inhibit  the  CuZnSOD  enzyme  activity. 

Dihydroethidium/dichlorofluorescein  diacetate  staining 

Superoxide-specific  dihydroethidium  (DHE)  analysis  was  per¬ 
formed  as  previously  described  [46].  Briefly,  superoxide  production 
was  estimated  using  the  fluorescent  dye  DHE  from  Molecular  Probes 
(Eugene,  OR,  USA).  Cells  were  washed  once  with  PBS  and  labeled  in 
suspension  at  37  °C  for  45  min  in  PBS  (containing  5  mmol/L  pyruvate) 
with  DHE  (lOpmol/L;  in  1%  dimethyl  sulfoxide  (DMSO)).  Samples 
were  analyzed  using  an  LSR  flow  cytometer  (superoxide-specific 
excitation  405  nm,  emission  585  nm,  band-pass  filter;  nonspecific 
excitation  488  nm,  emission  585  nm,  band-pass  filter).  The  mean 
fluorescence  intensity  (MFI)  of  10,000  cells  was  analyzed  in  each 
sample  and  corrected  for  autofluorescence  from  unlabeled  cells. 
Dihydrodichlorofluorescein  diacetate  (DCFH-DA)  analysis  was  per¬ 
formed  as  previously  described  [47].  Briefly,  steady-state  levels  of  pro¬ 
oxidants  were  determined  using  the  peroxide  oxidation-sensitive 
(CDCFH2,  lOpg/ml)  and  oxidation-insensitive  (CDCF,  lOpg/ml) 
fluorescent  dyes  (dissolved  in  1%  DMSO)  obtained  from  Molecular 
Probes.  The  cells  were  washed  once  with  PBS  and  labeled  with  the 
fluorescent  dyes  for  15  min  at  37  °C  in  PBS.  The  cells  then  were  then 
resuspended  in  PBS  and  analyzed  using  an  LSR  flow  cytometer 
(excitation  488  nm,  emission  530  nm,  band-pass  filter).  The  MFI  of 
10,000  cells  was  analyzed  in  each  sample  and  corrected  for 
autofluorescence  from  unlabeled  cells. 

Aconitase  and  succinate  dehydrogenase  activities 

Aconitase  activity  assay  was  adapted  from  the  previously 
described  experiment  [48].  Briefly,  protein  from  freshly  isolated 
thymocytes  was  isolated  by  sequential  freeze/thaws  (3x)  in  a  Tris/ 
MnCl2/citrate  buffer  (50  mM  Tris,  600  pM  MnCl2,  and  5  mM  Na 
citrate).  Protein  was  then  quantified  by  the  Bradford  assay.  After  this, 
200  pg  protein  was  combined  with  NADP+  (200  pM)  and  isocitrate 
dehydrogenase  (10  units).  Reaction  was  monitored  at  340  nm  for  the 
appearance  of  NADPH  formed  every  5  min  for  2  h.  Rates  were 
determined  from  slopes  determined  by  regression  analysis  of  data. 
Succinate  dehydrogenase  activity  assay  was  adapted  from  the 
previously  described  experiment  [49].  Briefly,  protein  from  freshly 
isolated  thymocytes  was  isolated  by  sequential  freeze/thaws  (3x)  in  a 
phosphate  buffer  (20  mM).  Protein  was  then  quantified  by  the 
Bradford  assay.  After  this,  200  pg  protein  was  combined  with 
succinate,  non-complex  11  inhibitors  (antimycin  A,  rotenone,  and 
cyanide),  CoQ,  and  2,6-dichloroindophenol  (DC1P).  Reaction  was 
monitored  at  600  nm  for  the  colorimetric  change  in  DCIP  accepting 
electrons  every  5  min  for  1  h.  Rates  were  determined  from  slopes 
determined  by  regression  analysis  of  data. 

Electron  microscopy 

Isolated  thymocytes  were  fixed  overnight  with  2.5%  glutaraldehyde 
in  0.1  M  cacodylate  buffer.  Postfixation  was  carried  out  for  1  h  at  room 
temperature  with  a  buffered  1%  osmium  tetroxide  solution  reduced 
with  1.5%  potassium  ferrocyanide.  Samples  were  stained  en  bloc  with 
2.5%  uranyl  acetate.  Cells  were  then  rinsed  and  dehydrated  using 


gradually  increasing  concentrations  of  acetone  to  100%.  Infiltration  of 
Spurr’s  epoxy  resin  and  acetone  was  carded  out  over  several  days  to 
100%  resin  and  cured  48  h  in  a  60  °C  oven.  Sections  of  90  nm  thickness 
were  cut  using  an  Ultracut  E  ultramicrotome  (Reichert-Jung).  Grids 
were  then  counterstained  with  5%  uranyl  acetate  for  2  min  and 
Reynold's  lead  citrate  for  2  min.  Samples  were  imaged  using  a  JEOL 
1230  transmission  electron  microscope  at  120  kV.  To  perform  quanti¬ 
tative  morphometry,  100  random  images  of  both  SOD2L/L  and  SOD2_/_ 
thymocytes  were  chosen  and  blindly  analyzed  using  ImageJ  (National 
Institutes  of  Health)  for  the  following  parameters:  cellular  area,  nuclear 
area,  cytoplasmic  area,  number  of  mitochondria,  mitochondrial  perim¬ 
eter,  mitochondrial  area,  and  mitochondrial  density. 

Flow  cytometry 

T-cell-specific  marker  analysis  for  development  or  post-influenza 
infection  was  performed  as  previously  described  [50],  After  isolation  of 
fresh  cell  preparations,  cells  were  suspended  in  staining  buffer  (HBSS 
supplemented  with  5%  bovine  calf  serum  and  0.1%  NaN3)  and  incubated 
with  conjugated  monoclonal  antibodies  in  the  presence  of  normal  rat 
serum  (to  limit  nonspecific  antibody  binding).  After  incubation  and 
washing,  the  cells  were  suspended  in  fixative  (1%  formaldehyde  in 
1.25x  PBS).  Stained  cells  were  run  on  a  FACSVantage  SE  flow  cytometer 
(Becton-Dickinson,  Mountain  View,  CA,  USA)  with  a  minimum  of 
30,000  events  collected  per  sample.  Antibodies  were  semipurified  from 
HB1 01  serum-free  supernatants  by  50%  ammonium  sulfate  precipitation 
and  conjugated  to  specific  fluorescent  dye  by  using  standard  procedures. 

Apoptosis  and  cell  cycle 

For  apoptosis,  the  Annexin  V-FITC  Apoptosis  Detection  Kit 
(Becton-Dickinson)  was  used  to  assess  annexin  V-positive  cells. 
Briefly,  fresh  cell  preparations  were  incubated  with  lx  annexin 
binding  buffer  and  annexin  V-FITC  (2.5  pg/ml)-conjugated  primary 
antibody  for  15  min  on  ice.  After  incubation,  propidium  iodide  (PI; 
10  pg/ml)  was  added  to  the  suspension  and  the  cells  were  analyzed  by 
flow  cytometry  using  an  LSR  flow  cytometer.  For  cell  cycle,  analysis 
was  adapted  from  a  previously  described  report  [51  ].  One  million  cells 
were  centrifuged  and  suspended  in  0.5  ml  of  Krishan  reagent  (0.1%  Na 
citrate,  0.03%  NP-40,  0.05  mg/ml  PI,  0.02  mg/ml  RNase  A)  before 
analysis.  Analysis  was  performed  by  examining  propidium  iodide- 
stained  nuclei  on  an  LSR  flow  cytometer. 

Influenza  and  rescue 

Mice  (weighing  approximately  20-25  g)  were  infected  with  1600 
TCIU50  of  mouse-adapted  A/PR/8/34  influenza  virus  A  H1N1  by 
intranasal  administration.  For  morbidity/mortality  experiments,  the 
mice  were  observed  and  weighed  daily  for  14  days.  For  pulmonary 
CD8+  T  cell  analysis,  mice  were  sacrificed  at  day  8  postinfection  and 
T  cells  isolated  from  lung  tissue.  1AV  PA224  (H2D(b)/SSLENFRAYV) 
and  NP366  (H2D(b)/ASNENMETM)  tetramers  were  obtained  from  the 
National  Institute  of  Allergy  and  Infectious  Disease  MHC  Tetramer  Core 
Facility  (Atlanta,  GA,  USA).  For  rescue  experiments,  Tempol  was 
obtained  from  Sigma-Aldrich  Co.  Mito-CTPO  was  graciously  donated 
by  Dr.  Balaraman  Kalyanaraman  (Medical  College  of  Wisconsin).  Both 
Tempol  and  Mito-CTPO  were  dissolved  in  ethanol  vehicle  and 
administered  at  25  or  250  nM,  respectively.  Water  bottles  were 
supplemented  with  sucrose  (4g/100ml)  to  offset  the  taste  of  the 
antioxidants.  Mice  were  treated  from  the  time  of  weaning  until  the 
date  of  experimentation. 

Statistics 

Data  are  expressed  as  means  and  standard  deviation.  All  experi¬ 
ments  were  performed  on  at  least  three  mice.  For  most  experiments, 
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comparisons  between  groups  were  analyzed  by  the  unpaired  two- 
tailed  Student  t  test.  For  Kaplan-Meier  analysis,  log-rank  analysis  was 
performed.  A  p  value  of  less  than  0.01  was  considered  significant. 


Results 

The  conditional  loss  of  SOD2  increases  superoxide-specific 
oxidative  stress  in  T  cells 

Because  of  the  lack  of  a  sufficient  model  to  study  ROS  stress  on  in  vivo 
development  of  the  mammalian  adaptive  immune  system,  we  reasoned 
that  the  creation  of  a  T-cell-specific  SOD2  knockout  mouse  would  serve 
as  an  optimal  platform  to  study  the  effects  of  mitochondrial  superoxide 
on  immune  system  development  and  function.  Floxed  SOD2  mice  were 
bred  with  mice  expressing  Cre  recombinase  under  the  control  of  the 
proximal  promoter  of  Lck  to  the  F2  generation  to  create  homozygous  T 
cell  SOD2  knockouts.  PCR  analysis  of  tail  sample  DNA  was  used  to 
confirm  the  proper  genotype  of  animals  used  for  studies  (Supplemental 
Figs.  1A  and  B).  After  this,  thymocytes  were  removed  from  animals  to 
confirm  the  recombination  of  SOD2  in  a  tissue-specific  manner. 
Genomic  DNA  extracted  from  T  cells  showed  the  specific  excision 
product  in  SOD2_/“  thymocytes,  with  no  apparent  Cre-recombination 
events  in  SOD2L/L  animals  (Fig.  1A).  Furthermore,  SOD2-specific  mRNA 
transcript  levels  were  significantly  decreased  in  knockout  animals 
compared  to  control  animals  (Fig.  IB).  T-cell-specific  SOD2~/_  mice 
were  found  to  have  minimal  SOD2  protein  and  activity  in  thymus- 
derived  T  cells  (Figs.  1C  and  D),  whereas  no  significant  decreases  were 
noted  in  any  off-target  tissues  (e.g.,  heart,  liver,  brain,  pancreas,  or 
muscle)  in  either  SOD2_/“  or  SOD2L/L  mice  (data  not  shown).  At  every 
level  of  analysis,  a  small  amount  of  residual  SOD2  was  detectable, 
suggesting  the  Cre  recombination  may  have  been  incompletely 
penetrant  or,  alternatively,  that  our  thymocyte  preparations  may 
contain  non-Lck-expressing  cells  (i.e.,  yd  T  cells,  immature  T  cell 
populations,  or  contaminating  stromal  cells)  and  thus  do  not  possess  Cre 
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Fig.  1.  SOD2  is  conditionally  and  effectively  removed  from  T  cells.  (A)  Conventional  PCR 
analysis  of  floxed  or  excised  SOD2  products  from  genomic  DNA  extracted  from  6-week- 
old  SOD21/I  orSOD2  1  thymocytes.  (B)  Quantitative  real-time  RT-PCR  analysis  ofSOD2 
mRNA  extracted  from  6-week-old  SOD2L/L  or  SOD2  thymocytes.  Results  are  shown 
as  normalized  to  SOD2L/L.  (C)  Comparison  of  SOD2  protein  in  6-week-old  SOD2L/L  or 
SOD2  7  thymocytes.  (J-Actin  (ACT(J)  is  shown  as  a  control  for  loading  and  transfer.  (D) 
Activity  assay  performed  on  protein  extracted  from  6-week-old  SOD2L/L  or  SOD2  1 
thymocytes.  At  least  three  mice  per  experiment  were  analyzed;  data  are  shown  as 
means  and  SD.  *p<0.01  by  Student's  t  test  versus  SOD2L/L. 


recombinase.  In  either  case,  taken  together,  these  data  support  the  T  cell 
specificity  and  efficiency  of  the  model  of  SOD2  knockout. 

To  confirm  that  SOD2“/_  mice  did  in  fact  harbor  increased 
oxidative  stress  compared  to  their  SOD2L/L  littermates,  we  performed 
DHE  stain  flow  cytometry  as  a  relative  measure  of  superoxide  content 
within  the  T  cell  populations  of  these  animals.  It  was  observed  that 
SOD2_/_  mice  demonstrated  significantly  increased  steady-state 
levels  of  superoxide  compared  to  control  SOD2L/L  mice  (Fig.  2A). 
Importantly,  no  change  was  noted  in  nonspecific  DF1E  oxidation  or  in 
the  oxidation  of  a  peroxide-sensitive  probe,  DCFH-DA  (Fig.  2B).  We 
confirmed  and  extended  these  findings  in  two  separate  enzyme 
activity  assays  that  are  highly  specific  to  disruption  by  superoxide 
(aconitase  and  succinate  dehydrogenase),  as  both  were  shown  to  be 
significantly  decreased  within  the  SOD2-/_  T  cells  with  no  significant 
change  in  protein  levels  (Figs.  2C  and  D).  It  was  previously  shown  that 
reactive  nitrogen  species  (RNS)  may  also  inactivate  these  enzymes 
[52],  but  probing  cellular  proteins  for  damage  by  RNS  showed 
nondetectable  levels  of  nitrotyrosine  (data  not  shown).  Overall,  it 
appeared  that  the  conditional  loss  of  SOD2  caused  a  specific  and 
significant  increase  in  steady-state  levels  of  superoxide  within  the  T 
cell  population  ofSOD2~/_  mice;  therefore  these  mice  should  serve  as 
an  effective  model  system  for  studying  the  effects  of  excess 
mitochondrial  superoxide  on  the  development  and  function  of  the 
mammalian  T  cell  adaptive  immune  system. 

Increased  mitochondrial  superoxide  levels  perturb  the  normal  development 
of  murine  T  cells 


To  understand  the  consequences  of  increased  steady-state  levels 
of  mitochondrial  superoxide  during  T  cell  growth  and  selection, 
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Fig.  2.  The  loss  of  SOD2  causes  superoxide-specific  oxidative  stress  in  T  cells.  (A) 
Quantification  of  total  (entire  bar)  and  superoxide-specific  (black  section  of  bar)  DHE 
staining  in  6-week-old  thymocytes.  (B)  The  peroxide-sensitive  probe  DCFH-DA 
demonstrated  no  significant  change  between  SOD2L/L  and  SOD2_/_  animals.  (C) 
Quantification  of  total  cellular  aconitase  and  succinate  dehydrogenase  (SDH)  activity  in 
6-week-old  thymocytes.  Black  bars  indicate  SOD2L/L  and  open  bars  indicate  SOD2_/_. 
(D)  Western  blots  of  the  iron-sulfur-containing  proteins  aconitase  2  (AC02)  and 
succinate  dehydrogenase  B  (SDHB)  demonstrated  no  significant  change  in  protein 
levels.  For  Western  blots,  p-actin  (ACT(3)  is  displayed  as  a  control  for  loading  and 
transfer.  At  least  three  mice  per  experiment  were  analyzed;  data  are  shown  as  means 
and  SD.  *p<0.01  by  Student's  t  test  versus  SOD2L/L. 
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quantitative  morphometry  was  performed  on  transmission  electron 
microscopy  (TEM)  images  of  developing  thymocytes  (Fig.  3A).  It 
was  observed  that  cellular  cross-sectional  area  was  significantly 
increased  within  SOD2“/_  T  cells,  but  no  change  was  noted  in 
nuclear  area  or  mitochondrial  size  (Fig.  3B;  Supplemental  Figs.  2A 
and  B).  Conversely,  SOD2-/~  thymocytes  demonstrated  an  approx¬ 
imate  doubling  of  the  number  of  mitochondria  per  cell,  but  these 
mitochondria  were  less  electron  dense  (Fig.  3B),  suggesting 
dysfunctional  or  degrading  organelles,  as  has  been  observed 
previously  in  oxidatively  stressed  systems  [53],  These  three 
mitochondrial  parameters,  decreased  electron  density,  increased 
number  per  cell,  and  altered  metabolism  (Fig.  2C),  all  further 
support  that  the  mitochondria  within  SOD2~/_  T  cells  are  defective. 


To  expand  upon  these  observations,  the  role  of  mitochondrial 
autophagy  was  explored.  It  was  shown  that  autophagy  was  in  fact 
increased  in  the  SOD2-/_  T  cells  as  shown  by  increases  in 
autophagic  vesicles  and  the  cleaved  form  of  LC-3  protein  (Figs.  3C 
and  D).  Overall,  the  loss  of  SOD2  seems  to  have  profound  effects 
upon  T  cell  ultrastructure  as  well  as  metabolic  function. 

In  addition  to  the  morphological  changes,  decreases  in  thymocyte 
numbers  as  well  as  peripheral  T  cells  were  noted  (Fig.  4A).  When 
total  thymocytes  were  quantified  at  various  ages  throughout  the 
mouse  immune  system  development,  an  approximately  twofold 
decrease  was  noted  in  SOD2_/_  mice.  In  contrast  to  this,  mice 
harboring  a  SOD2  knockout  from  an  earlier  stage  in  development 
(i.e.,  Vav-iCre-mediated  recombination)  demonstrated  no  significant 


Fig.  3.  Subcellular  and  mitochondrial  alterations  are  consequences  of  T-cell-specific  loss  of  SOD2.  (A)  Representative  TEM  images  of  SOD2L/1  (left)  and  SOD2  1  (right)  6-week-old 
thymocytes.  Open  arrowheads  indicate  mitochondria.  (B)  Quantitative  morphometry  of  thymocyte  cellular  area  (left),  mitochondrial  number  (center),  and  mitochondrial  density 
(right).  Each  "X”  represents  one  data  analysis  point;  open  bars  represent  means.  (C)  Representative  TEM  images  of  two  SOD2  1  6-week-old  thymocytes.  Solid  arrowheads  indicate 
apparent  autophagic  vesicles  commonly  seen  in  knockout  animal  thymocytes.  (D)  Western  blot  analysis  comparing  the  increase  in  cleavage  products  (i.e.,  LC3B-11)  of  the  autophagic 
marker  LC3.  (i-Actin  (ACT(i)  is  shown  as  a  control  for  loading  and  transfer.  At  least  three  mice  per  experiment  were  analyzed.  *p<0.01  by  Student’s  t  test  versus  SOD2L/L. 
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Fig.  4.  Significant!  cell  developmental  aberrations  are  observed  in  S0D2  animals.  (A)  Left:  total  thymic  cellularity  counts  at  various  ages  of  T  cell  development  for  SOD2L/L  (blue  circles) 
and  SOD2  1  (red  diamonds)  animals.  Open  bars  show  SOD2L/L  averages,  and  solid  bars  show  SOD  2  '  averages.  Middle:  total  cell  counts  for  three  peripheral  lymphoid  organs  at  6  weeks 
of  age.  Right:  total  thymic  cellularity  counts  for  6- week-old  mice  that  were  created  using  the  Vav-iCre  promoter.  No  significant  change  in  thymic  numbers  is  noted  when  knocking  out  SOD2 
at  a  different  stage  inT  cell  development.  (B)  Left:  representative  flow  scatter  diagrams  of  SOD2L/Land  SOD2  7  6- week-old  thymocytes  for  the  developmental  markers  CD4  and  CD8.  Right: 
quantification  ofCD4+  and  CD8+  6-week-old  thymocytes.  Colored  arrowheads  (left)  correlate  to  respective  colored  quantification  (right);  solid  bars  and  arrowheads  indicate  SOD2L/Land 
open  bars  and  arrowheads  indicate  50D2  (C)  Quantification  of  flow  cytometric  analysis  of  yli 1  T  cells  (left),  as  well  as  for  markers  indicative  ofT  cell  activation,  MEL14and  CD44,  on 

splenic  CD4+  cells  (middle)  and  CD8+  cells  (right).  In  all  graphs,  solid  bars  indicate  SOD2L/L  and  open  bars  indicate  SOD2  At  least  three  mice  per  experiment  were  analyzed;  data  are 
shown  as  means  and  SD.  *p<0.01  by  Student's  t  test  versus  SOD2L/L. 


decrease  in  thymic  cellularity  (Fig.  4A,  Supplemental  Fig.  3A).  This 
finding  suggests  that  various  stages  in  T  cell  development  may  be 
differentially  sensitive  to  the  presence  of  excess  superoxide. 
Moreover,  in  contrast  to  the  multiorgan  failure  noted  in  the 
aforementioned  constitutive  SOD2  knockout  mouse  [40],  our 
preliminary  studies  have  revealed  no  overt  phenotype  in  unstressed 
animals  harboring  either  liver-specific  or  mammary-specific  knock¬ 
outs  of  SOD2  (Supplemental  Figs.  3A,  B,  and  C).  Taken  together,  the 
identification  of  a  cellular  defect  in  the  T  cell  population  due  to  the 
loss  of  SOD2  at  a  later  (Lck-Cre-mediated)  rather  than  an  earlier 
(Vav-iCre-mediated)  stage  of  development  suggests  that  superoxide- 
mediated  effects  are  tissue  and  time  dependent  and  not  a  systemic 
non-specific  phenomenon.  Furthermore,  the  decrease  in  thymic 
cellularity  of  the  Lck-Cre  SOD2-/“  mice  was  exacerbated  by  the  fact 
that  SOD2_/_  mice  also  displayed  reduced  percentages  of  mature 
CD4+  and  CD8+  thymocytes  (Fig.  4B),  a  trend  that  was  also  observed 
in  peripheral  lymphoid  organs  as  well  (Supplemental  Fig.  4B).  In 


contrast,  the  most  immature  T  cell  populations  in  the  thymus  (i.e., 
CD4~,  CD8~),  yb  T  cells,  as  well  as  nonlymphocyte  populations,  were 
shown  to  be  static  or  increased  within  SOD2_/_  mice  (Figs.  4B  and  C, 
Supplemental  Fig.  4A).  This  is  most  likely  attributable  to  the  fact  that 
none  of  these  populations  activates  the  Lck  promoter  and  thus  they 
do  not  express  the  exogenous  Cre  recombinase.  Conversely,  these 
populations  may  be  exhibiting  compensatory  up-regulation  due  to 
the  decrease  in  mature  a[3  T  cells,  a  phenomenon  that  has  been 
described  in  other  models  ofT  cell  deficiencies  [42,54],  Interestingly, 
while  examining  peripheral  T  cells,  we  observed  that  SOD2_/_ 
lymphocytes  expressed  a  greater  proportion  of  Mell4“/CD44+  cells, 
suggesting  a  more  activated  state  among  otherwise  naive  T  cells 
(Fig.  4C).  This  finding  was  further  supported  by  microarray  analysis  that 
demonstrated  up-regulation  of  numerous  genes  involved  in  T  cell 
activation  in  the  thymocytes  of  T  cell  SOD2-/_  mice  (Supplemental 
Figs.  5A  and  B),  and  a  similar  effect  on  gene  regulation  by  mitochondrial 
superoxide  has  been  recently  reported  [55],  An  alternative  hypothesis  to 
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explain  this  apparently  activated  T  cell  phenotype  could  be  the  well- 
described  phenomenon  known  as  homeostatic  expansion,  a  compen¬ 
satory  mechanism  to  replenish  the  global  loss  of  T  cells  throughout  the 
SOD2_/~  mouse  [56].  Taken  together,  these  data  suggest  that  increased 
mitochondrial-derived  superoxide  alters  intracellular  signaling  within 
the  developing  thymocytes  and  may  inappropriately  mimic  the 
activating  conditions  of  an  infection,  thus  inhibiting  proper  T  cell 
function  during  true  infection. 

The  decrease  in  SOD2~,~  T  cell  numbers  is  due  to  increased  apoptosis 

During  T  cell  development  apoptosis  is  the  major  selective 
mechanism  for  lymphocyte  death  [57],  With  the  understanding 
that  apoptosis  may  be  initiated  through  cytochrome  c  release  from 
damaged  mitochondria  or  through  the  non-caspase-mediated  p53 
pathway  [58],  we  postulated  that  increased  apoptosis  due  to 
dysfunctional  mitochondria  may  be  the  mechanism  for  decreased 
thymocytes  in  the  SOD2_/_  mice.  When  flow  cytometry  for 
annexin  V/propidium  iodide  was  performed,  a  twofold  increase  in 
the  apoptotic  fraction  of  thymocytes  was  found  in  the  SOD2_/_ 
mice  (Fig.  5A).  In  addition,  an  approximate  twofold  increase  in 
thymocyte  sensitivity  to  oxygen  toxicity  in  vitro  as  well  as 
increased  immunoreactive  p53  was  also  noted,  further  suggesting 
increased  susceptibility  of  SOD2_/_  thymocytes  to  death  by 
apoptosis  (Supplemental  Figs.  6B  and  C).  In  addition,  propidium 
iodide  cell  cycle  analysis  demonstrated  a  greater  cycling  population 


in  SOD2~/_  mice,  which  supports  a  role  for  compensatory 
proliferation  in  more  primitive  CD3-/CD4-/CD8“  cells  as  previ¬ 
ously  discussed  (Supplemental  Fig.  6A).  When  mice  were  treated 
with  small-molecule  superoxide  scavengers  (i.e„  Tempol  [59]  or 
Mito-CTPO  [60])  added  to  the  water  supply  from  the  time  of 
weaning,  the  antioxidant  supplementation  led  not  only  to  the 
complete  rescue  of  the  SOD2_/~  phenotype,  but  also  to  increased 
thymocyte  numbers  in  control  SOD2L/L  animals  by  approximately 
twofold,  while  decreasing  their  apoptotic  fraction  by  50%  (Fig.  5B). 
There  was  a  strong  correlation  between  apoptotic  fractions  and 
total  cell  counts,  suggesting  that  apoptosis  is  the  primary 
mechanism  for  increased  T  cell  death  in  SOD2~/_  animals 
(Fig.  5B).  These  data  suggest  a  tight  relationship  between  the 
regulation  of  steady-state  mitochondrial  superoxide  levels  within 
the  normal  developing  thymus  and  immune  system  aberrations. 

Mitochondrial  superoxide  mediates  susceptibility  to  influenza  A  H1N1 

The  lack  of  proper  T  cell  development  observed  in  the  SOD2~/_ 
mice  strongly  suggested  an  immunocompromised  state.  To  test  this 
hypothesis,  mice  were  challenged  with  a  sublethal  dose  of  IAV  FilNl, 
a  pathogen  that  causes  a  T-cell-mediated  immune  response,  and  their 
ability  to  mount  an  immune  response  was  observed.  Surprisingly, 
100%  of  SOD2~/_  animals  succumbed  to  the  infection  and  died, 
whereas  only  a  small  percentage  of  SOD2L/L  animals  exhibited 
mortality,  and  furthermore,  SOD2-/~  mice  demonstrated  increased 
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Annexin  V 


Fig.  5.  Apoptosis  explains  the  decrease  in  SOD2_/_  thymocytes,  which  may  be  rescued  by  superoxide  scavenger  supplementation.  (A)  Representative  annexin  V  and  propidium 
iodide  flow  scatter  diagrams  of  6-week-old  thymocytes  from  vehicle-treated  SOD2L/L  (left),  vehicle-treated  SOD2_/_  (middle),  and  Mito-CTPO-treated  SOD2_/_  (right)  mice.  (B) 
Left:  total  thymic  cellularity  with  and  without  pharmaceutical  antioxidant  supplementation.  Middle:  quantification  of  annexin  V  positive  and  propidium  iodide  negative  thymocytes 
with  and  without  pharmaceutical  antioxidant  supplementation.  The  colors  are  coordinated  to  the  arrowheads  in  (A).  Right:  regression  analysis  of  total  thymic  cellularity  versus 
apoptotic  fraction.  At  least  three  mice  per  experiment  were  analyzed;  data  are  shown  as  means  and  SD.  *p<0.01  by  Student's  t  test  versus  vehicle-treated  SOD2L/L.  'Fp<0.01  by 
Student's  t  test  versus  vehicle-treated  SOD2_/_. 
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morbidity  (i.e.,  weight  loss)  with  no  recovery  compared  to  SOD2L/L 
mice  (Fig.  6A).  Moreover,  treatment  of  SOD2-/_  animals  with  the 
small-molecule  superoxide  scavenger  Tempol  (Supplemental  Fig.  7) 
or  Mito-CTPO  (Fig.  6A)  decreased  this  weight  loss  and  rescued  their 
ability  to  survive  the  influenza  infection  with  delayed  or  no  mortality 
noted.  When  examining  immunological  parameters  of  the  infection, 
we  once  again  found  that  SOD2_/_  mice  manifest  their  immune 
defect  with  decreases  in  not  only  total  T  cells,  but  also  IAV-specific 
CD8+  IAV-peptide-MHC-tetramer+  cells,  interferon-7  (IFN-y)+  cells, 
tumor  necrosis  factor  a  (TNF-a)+  cells,  and  dual-expressing  lFN-y+ 
TNF-ct+  cells  (Fig.  6B;  Supplemental  Fig.  8A).  Interestingly,  of  the 
SOD2_/_  T  cells  that  were  present  at  the  site  of  infection  the  quantity 
of  IFN-y  expression  per  cell  was  no  less  than  in  control  SOD2L/L  T  cells, 
and  TNF-a  expression  actually  increased  within  SOD2"/_  T-cells 
(Supplemental  Fig.  8B).  These  findings  are  consistent  with  other 
studies  of  immunodeficient  states  and  sepsis  in  which  T  cells 
compensate  by  overproducing  proinflammatory  markers  in  what  is 
known  as  a  “cytokine  storm”  [61]  and  may  in  part  explain  the 
apparently  activated  state  of  the  naive  T  cells  in  the  SOD2~/_  mice. 
Overall,  alterations  in  mitochondrial  superoxide  steady-state  levels 
seem  to  significantly  affect  the  outcomes  associated  the  mammalian 
adaptive  response  to  viral  infection. 

Discussion 

The  current  findings  suggest  an  excess  of  mitochondrial  superox¬ 
ide  at  the  stage  of  Lck  activation  of  development  tilts  the  intracellular 
redox  potential  toward  apoptosis,  resulting  in  the  decreased  thymic 


cellularity  observed.  A  relative  increase  in  the  normal  steady-state 
levels  of  superoxide  had  severe  consequences  on  the  immunocom- 
petence  of  the  mammalian  adaptive  immune  system.  This  novel 
finding  complements  the  well-established  role  of  ROS  in  the  innate 
immune  system  and  adds  yet  another  effect  of  superoxide  and  its 
derived  reactive  oxygen  intermediates  on  normal  tissue  growth  and 
development.  However,  further  studies  are  required  to  examine 
specific  signaling  pathways  that  are  affected  by  increased  mitochon¬ 
drial  superoxide.  It  would  be  predicted  that  redox-sensitive  signaling 
cascades  such  as  NF-kB  [17],  AP-1  [15],  or  JAK/STAT  [18]  may  be 
affected  and,  as  such,  may  in  part  in  explaining  the  pseudo-activated 
state  and  microarray  results  observed  in  this  study.  Furthermore,  the 
role  of  posttranslational  modifications  on  proteins  such  as  the 
hypoxia-inducible  factor  could  also  provide  valuable  information  on 
how  ROS  regulate  cellular  signaling  [62].  Moreover,  how  ROS  affect 
peripheral  lymphocyte  mitosis,  proliferation,  and  expansion  in 
antigen-unchallenged  and  -challenged  mice  also  warrants  further 
examination.  It  may  also  be  beneficial  to  understand  the  role  of 
compartmentalization  of  superoxide  (e.g.,  cytoplasmic,  extracellular), 
and  as  such  these  studies  should  be  followed  by  similar  investigations 
using  SOD1  [63]  or  SOD3  [64]  conditional  knockout  animals.  Because 
both  SOD1  and  SOD3  constitutive  knockout  animals  are  viable  and 
have  no  gross  pathologic  phenotype  compared  to  the  constitutive 
SOD2  knockout  mouse  [31,39,40,65],  it  would  be  postulated  that 
superoxide  compartmentalization  may  explain  this  phenomenon  at 
least  in  part.  We  observed  that  elevated  mitochondrial  superoxide 
levels  had  severe  effects  on  T  cell  development,  but  it  may  be 
postulated  that  elevated  superoxide  in  different  cellular  and  tissue 


Fig.  6.  SOD2  /  mice  are  immunologically  susceptible  to  influenza  virus,  but  may  be  rescued  by  superoxide-scavenging  pharmaceuticals.  (A)  Left:  Kaplan-Meier  analysis  of  mice 
succumbing  to  IAV  HI  N1  infection  during  a  2-week  period  after  viral  administration  with  and  without  Mito-CTPO  administration.  Right:  relative  weight  loss  of  mice  over  a  two-week 
period  after  infection  with  IAV  with  and  without  Mito-CTPO  supplementation.  Squares,  vehicle-treated  SOD2L/L;  triangles,  vehicle-treated  SOD2  '  ;  circles,  Mito-CTPO-treated 
SOD2  .  (B)  Left:  total  lymphocyte  counts  harvested  from  lungs  of  both  SOD2L/L  and  SOD2  '  animals  8  days  post-influenza  virus  infection.  Also,  flow  cytometric  analysis  of  CD8+ 
influenza-specific  NP366  and  PA224  MHC-tetramer-positive  cells  (middle)  and  analysis  of  NP366-  and  PA224-specific  IFN-7  /TNF-a:  (right)  CD8+  T  cells  isolated  from  lungs 
8  days  post-IAV  H1N1  infection.  At  least  six  mice  per  experiment  were  analyzed;  data  are  shown  as  means  and  SD.  For  weights  and  bar  graphs,  *p<0.01  by  Student's  t  test  versus 
SOD2l/l.  For  mortality,  *p<0.01  by  log-rank  analysis  versus  SOD2L/L. 
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compartments  may  lead  to  a  spectrum  of  developmental,  functional, 
or  disease  (e.g.,  cancer)  predispositions.  Last,  we  have  examined 
mitochondrial  superoxide  and  the  consequences  of  its  excess  on  the 
development  of  the  T  cell  adaptive  immune  system.  An  in-depth 
analysis  of  other  reactive  oxygen  species  (e.g.,  hydrogen  peroxide, 
peroxynitrite,  hydroxyl  radical)  may  prove  highly  informative  to 
further  understand  the  part  each  plays  in  the  development  and 
function  of  organ  systems.  These  data  could  serve  as  a  platform  for 
targeting  specific  immunodeficiencies  that  are  currently  not  well 
understood  and  may  be  due  to  alterations  in  redox  status  and  tailoring 
antioxidant  therapies  accordingly. 

In  this  study,  it  was  demonstrated  that  the  loss  of  immunocom- 
petence  in  a  SOD2~/_  background  could  be  rescued  using  a 
mitochondria-targeted  superoxide  scavenger,  an  example  of  the 
aforementioned  tailoring  of  specific  therapies  to  immunodeficiencies, 
Serendipitously,  it  was  shown  that  SOD2L/L  mice  with  fully  developed 
and  functional  adaptive  immune  systems  also  responded  to  the 
pharmaceutical  intervention.  These  mice  were  shown  to  have 
significantly  increased  numbers  of  thymocytes  compared  to  untreated 
control  animals,  and  this  seemed  to  be  due  to  a  decreased  apoptotic 
fraction.  At  first  glance,  it  seems  that  this  increase  in  T  cells  could 
potentially  be  a  boost  to  immune  system  prowess,  but  SOD2L/L  mice 
treated  with  superoxide  scavengers  showed  no  significant  increase/ 
decrease  in  the  ability  to  recover  from  influenza  A,  H1N1,  infection 
(data  not  shown).  This  finding  is  most  probably  attributable  to  the 
hypothesis  that  the  increase  in  T  cells  is  due  to  an  inhibition  of  proper 
apoptotic  death  of  cells  destined  to  die  in  the  thymus  (i.e., 
nonfunctional  cells,  cells  that  are  self-recognizing,  etc.)  as  opposed 
to  the  expansion  of  pathogen-specific  immune  cells.  Although  in  the 
short  term  this  had  no  observable  consequence,  it  is  speculated  that 
the  preservation  of  cells  predestined  for  death  could  lead  to  the 
potential  for  autoimmunity.  Additional  studies  will  be  needed  to 
confirm  this  potential  role  of  superoxide  in  the  development  of 
autoimmune  diseases  in  both  mouse  and  human.  In  this  study,  equal 
numbers  of  male  and  female  mice  were  used  for  all  assays,  but 
anecdotally  it  was  observed  in  many  cases  that  females  demonstrated 
more  pronounced  outcomes  due  to  the  effects  of  altering  steady-state 
superoxide  levels.  With  the  understanding  that  females  are  more 
prone  to  autoimmune  diseases,  as  well  as  the  recent  observation  of 
different  antioxidant  capacities  in  females  and  males  [66],  this 
preliminary  observation  may  warrant  further  examination  of  this 
phenomenon  in  similar  models  of  adaptive  immune  system  antiox¬ 
idant  alterations. 

In  conclusion,  we  have  used  a  previously  undescribed  model  of 
T-cell-specific  superoxide  scavenging  deficiency  to  examine  the 
development  of  the  adaptive  branch  of  the  mammalian  immune 
system  and  found  that  excess  mitochondrial  superoxide  late  in  T  cell 
development  is  detrimental  to  normal  T  cell  maturation.  In  contrast, 
early  developmental  loss  of  S0D2  activity  displayed  no  gross 
alteration  in  thymic  cellularity,  indicating  the  importance  of  the 
temporal  nature  of  the  loss  of  SOD2  function  on  thymic  developmen¬ 
tal  processes.  In  addition,  we  have  shown  that  the  loss  of  SOD2  does 
not  exert  a  similar  damaging  effect  in  development  or  maturation  of 
other  organs,  including  liver  and  mammary  glands,  supporting 
previous  findings  by  others  using  SOD2  floxed  models  [37,67],  At 
this  time  the  nature  of  the  differences  in  cell-type  specificity  and 
developmental  timing  of  superoxide  toxicity  are  poorly  understood, 
but  our  findings  warrant  further  investigation  into  possible  mecha¬ 
nisms  mediating  these  differential  sensitivities.  Using  the  appropriate 
Cre-recombinase  mouse,  this  model  could  be  expanded  to  examining 
other  elements  of  the  immune  system  (e.g.,  B  cells  [68],  neutrophils 
[69],  macrophages  [69]),  different  time  points  in  development 
[41,42,70],  or  other  organ  systems  altogether  [53,71,72].  In  addition, 
because  of  the  simplicity  of  the  breeding  scheme  the  floxed  SOD2 
mouse  may  be  applied  in  combination  with  other  protein  deletions  or 
additions  (e.g.,  p53  null  [73],  Myc  overexpression  [74,75])  to  examine 


the  effect  of  superoxide  oxidative  stress  in  alternative  immune  model 
systems.  Overall,  this  new  model  and  these  observations  serve  as  an 
initial  step  in  understanding  the  role  of  ROS  in  development  and 
function  of  mammalian  organ  systems. 
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